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Abstract 


The  purpose  of  this  research  was  to  investigate  the 
char acter i st i cs  of  a  small  and  inexpensive  pressure  transducer 
made  from  piezoelectric  film.  The  pressure  transducer  was 
designed  to  be  small  and  surface  mountable  so  that  the 
pressure  distribution  could  be  measured  on  an  aerodynamic 
body.  A  charge  reset  technique  was  used  to  overcome  the 
traditional  1  imitations,  low  frequency  measurements  and  charge 
decay,  of  a  piezoelectric  device.  ^The  resulting  time  history 
had  to  be  integrated,  so  a  method  was  developed  to  integrate 
the  data  and  remove  sources  of  error.  The  time  history 
reconstruction  technique  was  tested  with  a  known  input  and 
found  to  accurately  reassemble  the  original  input.  The 
piezoelectric  film  pressure  transducer  was  tested  in  a  flow 
field  along  side  a  commercial  pressure  transducer .  The 
piezoel ectr  ic transducer's  frequency  response  showed  agreement 
with  the  frequency  response  of  the  commercial  transducer. 


Also,  the  time 

h  i  story 

response 

of 

the  piezoelectric 

transducer  gave 

ev i dence 

that 

1  ow 

frequency  pressure 

measurements  are 

poss i b 1 e  with  a 

piezoelectric  pressure 

transducer . 


V  7  /y  . 
/ 


x  i  i  i 


•  ■  W\  <* 


Development  of  a  Piezoelectric  Polymer  Film 
Pressure  Transducer  for  Low  Frequency  and 
Dynamic  Pressure  Measurement 


I .  I ntroduct i on 

Mot i vat i on 

With  increasing  cost  for  full-scale  testing,  more  model 
testing  and  simulation  is  being  done.  Some  of  the  typical 
wind  tunnel  data  collected  are  pressure  and  force 
measurements.  The  instrumentation  or  transducers  needed  to 
make  these  measurements  requires  input  and  output  through 
tubing  and  wires.  In  some  cases,  an  accurate  model  may 
require  many  transducers.  As  a  result,  the  transducers 
occupy  a  large  portion  of  the  model.  If  high  resolution 
measurements  are  needed,  a  model  sufficiently  large  to 
encompass  al  1  the  transducers  may  be  bigger  than  the  avai lable 
testing  facility.  To  accurately  obtain  data,  either  the  test 
facility  would  need  to  be  larger  or  the  model  would  need  to  be 
smaller  and  more  tests  run.  These  facts  increase  the  cost  and 
complexity  of  testing  the  model. 

If  a  smaller  and  cheaper  pressure  transducer  were 
available,  model  size  and  cost  could  be  reduced.  Potentially, 
a  small  surface  mounted  pressure  transducer  could  be  developed 
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to  cover  a  model  in  adequate  density  to  map  the  entire  surface 
pressure  field.  To  date,  no  transducer  has  been  developed 
that  is  small,  surface  mountable  and  inexpensive. 

Summary  of  Current  Knowledge 

In  the  late  1  960's  a  new  piezoelectric  polymer  film 
material  was  developed  (1:328).  Researchers  found  that 
po  1  y v i ny 1 i dene  fluoride  (PVDF)  could  be  processed  to  act  as  a 
piezoelectric  material.  Piezoelectric  materials  convert 
physical  motion  or  mechanical  stress  to  an  electrical  signal; 
the  inverse  process  can  also  be  realized  (1:328-329;  4:22-33). 
Since  pressure  exerts  a  force  on  the  body  on  which  it  acts, 
piezoelectric  film  could  be  used  as  a  pressure  transducer 
(5:5).  The  piezoelectric  film  would  output  an  electrical 
signal  in  response  to  a  pressure. 

Since  the  inception  of  PVDF  piezoelectric  film  in  the 
late  1960's  (  1:328-329),  it  has  been  used  in  hundreds  of 
different  applications.  Piezoelectric  film  has  been  used  to 
detect  heat,  act  as  a  switch,  to  pick  up  sound  from  musical 
instruments,  as  a  pressure  sensor,  and  even  to  generate 
electricity  (4:22-33).  As  a  pressure  sensor,  PVDF  film  has 
found  applications  by  organizations  such  as  NASA,  the  US  Navy, 
German  and  American  researchers,  and  many  others. 

In  the  mid  1970's,  the  National  Bureau  of  Standards  ( NBS ) 
developed  a  PVDF  transduce-"  to  study  impact  pressures  in  car 
crashes  (3).  They  developed  two  transducers:  1  )  a  1  cm 
diameter  sensor  and  2)  a  smaller  2  mm  diameter  sensor.  The 
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1  cm  sensor  worked,  but  it  produced  extraneous  outputs  due  to 
film  flexing  and  stretching  from  an  improper  mounting 
technique.  The  2  mm  diameter  sensor  provided  a  signal-to- 
noise  ratio  larger  than  the  1  cm  diameter  sensor.  The  stud> 
did  not  proceed  with  size  optimization,  but  recommended  that 
this  area  should  be  investigated. 

NASA  tested  PVDF  film  for  use  as  a  flow  diagnostic  tool 
(2).  They  investigated  two  sensor  configurations,  a  rigidly 
mounted  one  and  another  mounted  over  an  air  filled  cavity. 
Evaluation  of  the  sensors  included  sensitivity  and  frequency 
response  and  data  optimization.  The  results  demonstrated  that 
the  rigidly  mounted  piezoelectric  film  was  affected  by 
structural  responses;  whereas,  the  cavity  mounted  sensor  was 
influenced  significantly  less  by  the  structure.  Overall,  the 
cavity  mounted  sensor  exhibited  a  broader  and  more  sensitive 
frequency  response. 

German  researchers  used  p i ezoe 1 ectr i c  film  and  multi-sensor 
piezoelectric  arrays  to  measure  forces  in  flow  experiments 
(6).  They  developed  three  arrays  for  measuring  dynamic 
pressure  distributions.  The  arrays  were  shaped  as  strip, 
raster,  and  matrix  arrays.  The  researchers  conducted  wind 
tunnel  and  free  flight  tests  using  the  arrays  and  found  that 
the  array  concept  worked.  Only  basic  data  collection 
techniques  were  given;  further  development  would  be  required 
to  use  the  arrays  on  scale  models. 

Several  university  researchers  developed  a  new  charge 
readout  technique  for  low-frequency  pressure  measurements  with 
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piezoelectric  film  (7).  This  technique,  called  the 
incremental  Readout  Algorithm  (|RA),  operates  on  a  sample-and- 
reset  principle  that  overcomes  the  charge  decay  problem  of 
p  i  ezoe 1 ectr i c  devices.  The  IRA  allows  low-frequency 
measurements  with  a  high  sensitivity  and  dynamic  range.  These 
measurement  techniques  removed  the  low-frequency  limitations 
imposed  by  piezoelectric  film  pressure  transducers. 

Problem  Statement 

The  purpose  of  this  study  was  to  develop  a  small  and 
inexpensive  pressure  transducer  from  piezoelectric  film.  The 
transducer  would  be  capable  of  measuring  static  and  dynamic 
pressure.  In  addition,  algorithms  and  electrical  hardware 
were  developed  to  measure  and  reduce  the  data  collected  from 
the  transducer. 

Assumpt i ons 

The  following  assumptions  were  made: 

1.  Pyroelectric  effects  are  negligible  compared  to  the 
piezoelectric  response,  since  the  operating  temperatures  are 
stable  and  below  the  non-linear  temperature  region  of  the 
piezoelectric  film. 

2.  All  surface  stress  effects  are  small  and  can  be 
ignored  because  the  film  is  rigidly  attached  to  the  model’s 
surface . 


1-4 


Scope 


The  scope  of  this  study  was  to  design,  develop,  and 
evaluate  performance  of  a  piezoelectric  film  pressure 
transducer.  This  process  included  developing  digital  circuits 
to  collect  and  measure  the  transducer’s  response.  Algorithms 
were  developed  to  reduce  the  data  to  a  practical  (useable) 
format.  A  method  was  developed  to  apply  the  PVDF  film  to  a 
body  so  that  surface  pressure  could  be  measured.  The  research 
did  not  include  film  size  and  shape  optimization,  nor  were 
different  film  adhesives  addressed. 

The  study  used  commercially  available  components  and  did 
not  develop  any  new  items.  Using  existing  resources  helped 
keep  the  design  simple  and  inexpensive.  The  PVDF  film  was 
modified,  but  its  basic  physical  and  electrical  structure  was 
not  changed. 

The  project  concluded  when  a  pressure  measurement  was 
obtained  with  the  piezoelectric  transducers  mounted  on  a  body. 
Approach 

The  transducer  development  was  divided  into  three  phases 
which  included  the:  1)  development  of  a  charge  amplifier 
circuit  with  a  reset  capability;  2)  optimization  of  the 
charge  amplifier  to  minimize  DC  offset  and  noise;  and  3) 
performance  evaluation  of  the  transducer  in  a  dynamic 
situation  by  comparing  it  to  a  known  transducer.  Each  phase 
was  subdivided  into  incremental  steps,  they  are  listed  below: 
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1.  Phase  1  -  Development  of  a  charge  amplifier  with  reset 
capab i 1 i ty . 

o  Develop  the  electronics  required  to  discharge  the 
piezoelectric  film  charge, 
o  Develop  software  to  reduce  the  data. 

o  Refine  the  hardware,  software  and  procedures  to  obtain 
accurate  and  repeatable  data. 

2.  Phase  2  -  Optimization  of  the  charge  amplifier  circuit, 
o  Minimize  charge  leakage. 

o  Minimize  DC  offset  and  noise. 

o  Determine  the  upper  and  lower  reset  rate  limits, 
o  Characterize  the  charge  amplifier  response. 

3.  Phase  3  -  Development  of  a  film  application  for  pressure 
measurement . 

o  Determine  a  way  to  bond  the  PVDF  film  to  a  body, 
o  Test  the  film's  capability  to  measure  pressure. 

Order  of  Presentation 

Background  information  on  piezoelectricity  and 
piezoelectric  materials  are  presented  in  Chapter  II.  The 
third  chapter  presents  the  details  of  the  reset  scheme,  time 
history  r econstr uct i on  technique,  and  testing.  Chapter  IV 
discusses  test  and  evaluation  results.  The  last  chapter, 
Chapter  V,  summarizes  the  results  and  recommends  areas  of 
further  research  and  development  on  the  piezoelectric  film 
pressure  transducer. 
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Background 


History 

A  strange  phenomenon  was  long  known  by  the  people  in 
Ceylon  and  India.  When  a  tourmaline  crystal  was  thrown  into 
a  fire,  it  would  attract  the  hot  ashes  and  repel  them  when  it 
cooled.  In  1703,  Dutch  traders  brought  a  "Ceylon  Magnet"  back 
from  Ceylon  to  Europe  (8:1).  Linnaeus  gave  the  scientific 
name  Lapis  electricus  to  this  phenomenon  in  1747  (8:1).  In 
1756  Aepinus  investigated  this  effect  and  noted  the  electrical 
property  of  opposite  polarities  in  the  heated  tourmaline 
crystal  (8:1).  Another  investigator,  Brewster,  assigned  the 
term  "pyroe 1 ectr i c i ty "  to  this  phenomenon  after  he  had 
observed  it  in  tourmaline  and  other  crystals  (8:1).  The  term 
"pyro-"  comes  from  the  Greek  word  meaning  "fire",  and  the 
pyroelectric  effect  manifests  the  production  of  electricity  in 
a  material  caused  by  the  thermal  stresses  due  to  temperature 
changes . 

Other  investigators  continued  to  observe  the  pyroelectric 
effect  of  crystals.  But  in  1880,  the  Curie  brothers,  Pierre 
and  Jacques,  discovered  that  when  certain  crystals  were 
compressed  in  a  particular  direction,  a  charge  accumulation 
results  (8.2).  The  Curie  brothers  wrote  of  their  discovery  in 
a  paper  saying: 

Those  crystals  having  one  or  more  axes  whose 
ends  are  unlike,  that  is  to  say,  hemihedral 
crystals  with  oblique  faces,  have  the  special 
physical  property  of  giving  rise  to  two  electric 
poles  of  opposite  signs  at  the  extremities  of  these 
axes  when  they  are  subjected  to  a  change  in 
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temperature:  this  is  the  phenomenon  known  under 
the  name  of  pyroe 1 ectr i c i ty  .  We  have  found  a  new 
method  for  the  development  of  polar  electricity  in 
these  same  crystals,  consisting  in  subjecting  them 
to  variations  in  pressure  along  their  hemihedral 
axes .  (8:2) 

Fol lowing  their  initial  discovery,  the  Curies  observed  the 
"polar  electric"  effect  in  other  crystals  such  as  :  zinc 
blende,  sodium  chlorate,  boracite,  quartz,  calamine,  topaz, 
tartaric  acid,  cane  sugar,  and  Rochelle  salt  (8:3). 
Immediately,  a  lot  of  interest  was  aroused  in  the  scientific 
community.  One  of  the  interested  parties,  Hankel,  proposed 
the  name  "piezoelectricity",  since  he  thought  the  new 
discovery  followed  its  own  unique  laws  and  was  unrelated  to 
py roe  1 ectr i c i ty  (8:3).  From  the  Curies’  discovery,  arose  both 
the  name  and  field  of  piezoelectricity.  Similar  to  the  term 
"pyro-",  "piezo-"  means  "to  press"  in  the  Greek  language  and 
the  piezoelectric  effect  is  caused  by  applying  pressure  to  the 
material.  In  1881,  the  converse  of  piezoelectricity  was 
proposed  by  Lippman,  and  the  Curies'  verified  this  the  same 
year  (8:4). 

Many  others  began  to  investigate  the  piezo-  and  pyro¬ 
electric  phenomena.  Among  the  most  noted  was  Woldmar  Voight, 
who  based  his  theory  on  the  thermodynamic  principles  of  Lord 
Kelvin  (8:5).  In  1910,  Voight  published  a  work  called 
L ehrbuch  der  Kr i sta 1 1 phys i k ,  which  became  a  standard  reference 
for  p i ezoe  1  ectr i c  researchers  (8:5).  His  work  identified  32 
crystal  classes  that  could  display  p i ezoe 1 ectr i c  properties; 
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Figure  2-1.  Barium 
Titanate:  A  Piezoelectric 
Crystal  (9:1758). 


and  of  the  18  piezoelectric  coefficients,  he  showed  which  ones 
had  nonzero  values  (8:5). 

For  the  next  three  decades,  piezoelectricity  was  confined 
to  research,  with  no  practical  application.  In  France,  at  the 
beginning  of  World  War  I,  a  scientist  by  the  name  of  Longevin 
developed  a  means  to  piezoelectrical ly  exciting  quartz  plates. 
These  plates  could  be  used  as  high  frequency  emitters  and 
receivers  of  sound  waves  in  water.  Longevin  and  others  had 
developed  an  echo  method,  SONAR,  which  could  locate  objects 
underwater  and  map  the  contour  of  the  ocean  floor  (8:5). 
Thus,  Longevin  became  the  father  of  ultrasonics,  which  is 
widely  used  today. 

During  World  War  I,  research  discovered  many  uses  of 
piezoelectricity  and  the  properties  of  crystals.  Some  of  the 
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more  useful  purposes  were;  Rochelle  salt  crystals  used  for 
underwater  signaling  and  quartz  used  for  precision  cIocks 
(8:6).  Researchers  at  MIT,  in  the  late  1930's  found  that 
certain  ceramics  such  as  barium  titanate  (see  Figure  2-1)  and 
lead  zirconate  titanate  (PZT)  could  be  polarized  to  respond 
piezoelectrical 1y  (7:1).  Researchers  continued  to  investigate 
other  p i ezoe 1 ectr i c  materials  such  as  ceramics,  crystals, 
organic  materials,  human  tissue,  and  bone  (7:1). 

In  the  late  1960's,  a  Japanese  researcher,  Kawai,  found 
that  a  polarized  polymer,  po 1 y v i ny 1 i dene  difluoride  (PVDF), 
developed  a  much  greater  piezoelectric  charge  compared  to 
other  polymers  (7:1,  1:328).  This  began  the  vast  field  of 
development  of  piezoelectric  film. 

P i ezoe 1 ectr i c i tv 

Before  delving  into  a  discussion  of  the  fundamentals  of 

piezoelectricity,  two  definitions  by  W.G.  Cady  are  useful. 

Cady  defines  the  piezo-  and  pyro-  effects  as  follows: 

Piezoelectricity  is  electric  polarization  produced 
by  mechanical  strain  in  certain  crystals,  the 
polarization  being  proportional  to  the  amount  of 
strain  and  changing  sign  with  it.  The  reverse  is 
also  true;  that  i  s .  . . an  electrical  polarization 
will  induce  a  mechanical  strain  in  p i ezoe 1 ectr i c 
crystals. 

Pyroelectricity  is  electric  polarization  induced  by 
thermal  absorption  in  certain  crystals,  the 
polarization  being  proportional  to  the  level  of 
thermal  change  (8:4). 

Piezoelectric  Effect.  The  piezoelectric  effect  is  only 
found  in  materials  that  have  "ionic”  bo.  ds  (atoms  that  are 
arranged  in  positive  and  negative  pairs  called  dipoles) 


2-4 


(9:1757).  When  pressure  is  applied  to  the  material,  the 
dipoles  are  distorted  or  polarized.  This  behavior  causes  a 
separation  about  the  center  of  inversion  which  creates  an 
electric  field  (see  Figure  2-2).  The  center  of  inversion,  or 
geometric  center  gives  rise  to  the  property  of 
piezoelectricity.  If  a  crystalline  material  has  perfect 
symmetry,  no  combination  of  uniform  stresses  will  produce  the 
distortion  required  for  the  piezoelectric  effect.  A  common 
crystal  with  a  perfectly  symmetrical  structure  is  sodium 
chloride  (or  salt),  each  positive  charge  is  surrounded 
symmetrically  by  a  negative  charge.  Barium  titanate  (see 
Figure  2-1)  possesses  an  asymmetrically  crystalline  structure 
and  is  piezoelectric.  It  is  one  of  the  20  (of  the  32  overall) 
crystal  classes  that  possess  p i ezoe 1 ectr  i  c  properties  due  to 
the  asymmetries  (9:1757). 

Dipoles.  The  electric  dipole  helps  explain  why 
piezoelectricity  occurs  on  a  microscopic  basis.  An  electric 
dipole  consists  of  two  point  charges  of  equal  magnitude,  but 
different  sign,  which  are  spaced  at  a  fixed  distance  (10:132). 

Figure  2-3  shows  two  equal  and  opposite  fundamental  units 
of  charge  q  spaced  a  distance  d  apart  with  a  dipole  strength 
of  pq  (10.132).  The  electrical  potential  <t>  at  an 
arbitrary  point  p,  caused  by  +q  and  -q  are  given  in  equations 
2-1  and  2-2. 
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Figure  2-2.  Concept  of  a  Center  of  Inversion  in  a  Crystalline 
Structure.  (a)  Crystal  with  a  Center  of  Inversion,  and  (b) 
Crystal  without  a  Center  of  Inversion  (24:407). 
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Where  €,  is  the  permittivity  of  vacuum.  The  total  electric 
potential  <j>  at  p  is  the  superposition  of  the  two  potentials 
caused  by  the  individual  positive  and  negative  charges.  That 

1S  ’  /  \ 

♦  +  •  (2_3) 
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Figure  2-3.  Simple  Dipole  (16). 


If  r  is  large  compared  to  d,  the  distances  from  q  to  p  can 
be  simplified  with  the  following  approximations  (10.133). 
That  is, 

r1=r+(~)c°s0  (2-4) 

r2=r-|-^|cos0  .  (2-5) 

The  total  potential  at  point  p  becomes: 


q  d  cosfl  > 

4ue0ra  ..  I  d  cos8\2 
\  \  2  r  }  ) 


(2-6) 


Since  r  was  assumed  to  be  much  greater  than  d,  the  denominator 


can  be  approximated  with  the  following  form: 


(j  (j  CO30 

4ne0ra 


(2-7) 
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If  qd  is  defined  as  the  dipole  moment  P  (P  =  qd),  the 
numerator  of  equation  2-7  can  be  written  as  a  dot  product; 


that  is,  qdcosO  =  P»r  (10.133). 

The  electric  field  E  of  the  dipole  is  obtained  from  the 
relation  E  =  -V<j>.  Thus,  the  general  form  of  dipoles  can  be 
written  in  vector  form  as  shown  in  equation  2-8.  That  is, 


(2-8) 


To  consider  a  simple  example  in  a  crystal,  the  electric  field 
is  expressed  below  in  polar  coordinates  (16:36).  That  is, 

E- — — cos0  f  +  — sinB  0  .  (2-9) 

2  sear3  4  *e0r3  v  ' 


In  an  equilibrium  or  unstressed  state  (see  Figure  2-4a)  ,  there 
are  no  deflections  or  electric  fields.  All  dipoles  are  in  the 
same  plane,  equidistant,  and  symmetrical  about  the  negative 
charge,  -q.  The  total  electric  field,  Et#ui ,  is  the  sum  of  all 
contributions  to  the  field  (11:366).  That  is, 

£**•2=5 ’l+^+*3  •  (2-10) 

Since  the  dipoles  are  symmetrical  and  equal  in  magnitude  about 
-q,  all  the  dipoles  contribute  equally  and  cancel  their 
overall  effect.  Thus, 

and  EtotAl- 0  .  (2-11) 

When  a  stress  is  applied,  as  shown  in  Figure  2-4b,  the 
noncentral  ions  displace  angularly  about  the  center,  creating 
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1 20c 

1  30' 

(  o  ) 

(  b  ) 

Figure  2-4.  Example  of  (a)  Unstressed  and  (b)  Stressed  Crystal 
with  Angular  Rotations. 

an  electric  field.  In  this  case,  the  individual  contributions 
are  unequal.  That  is, 

Ex  *  E2  *  Ej  .  (2-12) 

The  resulting  expression  for  EtsU(  is  shown  below.  That  is, 

E=E1+E3+E3  =  -^=L  (0.155)0  (2-13) 

where  d  is  the  radial  distance  from  -q  to  the  positive 
charges,  and  r  is  the  distance  to  an  arbitrary  observation 
point. 

Though  the  microscopic  perspective  contributes  an  insight 
into  the  cause  of  the  p  i  ezoe  1  ectr  i  c  effect,  it  does  not 
adequately  describe  the  microscopic  behavior  of  crystals. 
Cady  states  the  limitations  of  this  atomic  theor y .  That  is. 
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In  spite  of  the  fact  that  molecular  or  atomic 
theories  of  p i ezoe 1 ectr i c i ty  began  to  appear  very 
soon  after  the  Curies'  discovery,  a  satisfactory 
theoretical  treatment  of  the  phenomenon  can  hardly 
be  said  to  have  passed  the  initial  stage.  The 
resources  of  modern  lattice  dynamics  are  still 
unequal  to  the  task  of  predicting  anything  better 
than  a  rough  approach  to  the  order  of  magnitude  of 
the  piezoelectric  effect,  even  for  the  simplest 
structures.  Incomplete  though  the  story  is, 
however,  it  offers  much  of  interest  on  both  the 
theoretical  and  the  experimental  side.  (8:731) 


Piezoelectric  Properties 

The  microscopic  perspective  of  the  p  i  ezoe 1 ectr i c  effect 
helps  to  explain  why  a  charge  develops  in  a  p i ezoe 1 ectr  i  c 
material,  but  it  does  not  explain  the  macroscopic  or  the 
observable  properties.  Generally,  a  book  would  need  to  be 
written  to  completely  describe  all  the  properties  and  their 
applications;  only  a  brief  overview  of  selected  ones  will  be 
given. 

Piezoelectric  materials  generate  a  charge  when  they  are 
deformed.  This  effect  is  reversible  because  if  a  charge  is 
applied  to  the  material,  it  will  respond  with  a  mechanical 
deformation  (12:266).  The  electro-mechanical  energy 
conversion  principle  is  very  useful,  and  it  is  the  basis  of 
many  instruments  that  measure  force,  pressure,  and 
acceleration  (12:267). 

D i r ect i ona 1 i ty .  Piezoelectric  materials  display 
anisotropic  behavior,  because  they  respond  differently, 
depending  upon  the  material 's  orientation.  The  electrical  and 
mechanical  responses  differ  according  to  the  axis  that  the 
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electric  field  or  mechanical  stress  or  strain  is  applied 
(5:9).  In  making  any  calculations,  the  material's 
directionality  must  be  accounted  for. 

The  sign  convention  used  throughout  this  development  1  s 
shown  in  Figure  2-5.  The  material's  length  1,  corresponds  to 
the  1-axis,  2  to  the  width  w,  and  3  to  the  thickness  t  (5:9). 
Double  subscripts  are  used  on  numerical  physical  constants 
(presented  in  the  following  sections)  to  describe  the 
occurring  phenomenon  (5:9;  12:267).  The  first  subscript 
identifies  the  axis  or  direction  of  the  electrical  effect. 
The  second  subscript  refers  to  the  axis  of  the  mechanical 
effect  (12:267). 


3 


Axis  System  of  a  Piezoelectric  Material 


Figure  2-5.  Axis  System  for  a  P i ezoe 1 ectr i c  Material  (5:9). 

Tensile  forces  or  strains  imply  positive  mechanical 
action,  and  compression  is  negative.  Positive  electrical 
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effects  result  in  an  increase  in  the  charge  polarization; 
negative  electrical  effects  cause  a  decrease  in  the  charge 
polarization  (5:9). 

Curie  Temperature.  The  Curie  temperature  is  the 
temperature  above  or  below  which  the  dipoles  that  compose 
electrets  become  randomly  oriented  (12:267).  As  a  result,  the 
material  looses  its  piezoelectric  property.  If  the  material's 
temperature  exceeds  the  Curie  point,  it  will  loose  its 
piezoelectric  effect;  upon  cooling,  it  returns  to  being 
piezoelectric  only  if  it  is  poled  again.  Thus,  the  useful 
range  of  a  piezoelectric  material  is  below  the  Curie  point 
(9:1758).  The  Curie  point  can  also  be  changed  in  some 
materials  by  changing  their  chemical  composition;  this  feature 
results  in  a  more  satisfactory  Curie  temperature  (9:1758). 

Capac i tance .  Typically,  piezoelectric  materials  are 
inherently  capacitive.  In  use,  metal  electrodes  are  deposited 
on  two  opposing  faces  of  the  material  and  wires  are  attached 
to  access  the  developed  charge.  Piezoelectric  materials  are 
insulators  and  the  electrodes  become  capacitor  plates. 
Mechanical  distortions  result  in  a  charge  which  in  turn 
appears  as  a  voltage  between  the  electrodes.  The  resulting 
voltage  follows  the  law  for  capacitors  shown  in  equation  2-14. 

E  -2  (2-14) 

where  E  is  the  voltage,  Q  is  the  charge,  and  C  is  the 
capac i tance . 
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Electromechanical  Coupling  Factor.  The  coupling  factor 
K*  defines  a  piezoelectric  material's  ability  to  exchange 
electrical  and  mechanical  energy,  and  the  reverse  (5:18).  For 
example,  k:,  would  describe  the  electrical  energy  brought 
about  by  a  mechanical  strain  divided  by  the  input  electrical 
energy  (5:18).  The  following  relations  describe  k! : 

£2  -  mechanical  energy  converted  to  electrical  energy 

input  mechanical  energy  (2-15a) 

£2  -  electrical  energy  converted  to  mechanical  energy 

input  electrical  energy  .(2-1 5b) 

Mechanical  Constants.  There  are  three  main  types  of 
mechanical  constants  for  p i ezoe 1 ectr  i  c  materials.  They  are 
the  strain  constant,  d,  the  stress  constant,  g,  and  a  related 
stress  constant,  e.  The  g  constant  gives  the  ratio  of  the 
electric  field  generated,  under  open  circuit  conditions,  to 
the  stress  that  is  applied  along  a  particular  axis.  The  same 
constant,  g,  also  gives  the  proportional  amount  of  strain 
induced  under  stress  free  conditions,  along  an  axis  compared 
to  the  electrical  charged  applied  per  unit  area  of  the 
electrode  (5:17).  Equations  2-16a  and  2-16b  show  g  i  r. 
equation  form,  where  D  is  the  electrical  displacement,  E  the 
electric  field,  X  the  stress,  and  x  the  strain. 

electrical  field  developed  _/  (2-16a) 

9  applied  mechanical  stress  \^xj  D  '  ’ 


strain  developed  _/  dx\ 

9  ~  appl  i ed  charge  densi  ty  \dD}x 


( 2- 1 6b ) 
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The  strain  constant,  d,  gives  the  change  in  the 
piezoelectric  material  dimensions,  in  all  three  axes,  when 
electric  field  is  applied  to  the  material  in  a  stress  free 
state  (5:17).  Thus,  d  gives  the  proportionality  constant  of 
the  stress  free  dimensional  change  relative  to  the  field 
applied.  Additionally,  d,  (see  Figure  2-6)  expresses  the 
ratio  of  the  short  circuit  charge  per  unit  area  of  the 
electrode  to  the  stress  applied  along  a  particular  direction 
(5:17).  These  relationships  are  expressed  in  the  following 
equat ions: 

d_  electrical  charge  denai  ty  developed  _  /  dD\  ,  ) 
applied  mechanical  stress  \dx)s  ' 


d= - strain  developed  /to 

appl i ed  el ectzi cal  field  \  dB) x 

The  constants  can  be  calculated  from  both  relationships 
provided  the  dielectric  constant,  e,  of  the  material  is  known. 
That  is, 


C- 


ewl 

t 


(2-18) 


where  w  is  the  width,  t  is  the  thickness,  1  is  the  length,  and 
C  is  the  capacitance.  Finally, 


electrical  field  developed  _  eow^  _  eoc  _  _Q_  _  d 
applied  mechanical  stress  tf  cf  ef  e 


(2-19) 


or 


d-eg 


(2-20) 
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where  e„  A  electrical  field  developed 


f  A  app 1 i ed  force 
Q  A  charge  deve 1  oped (  1  2 : 268 ) . 


x2=x3=0  l  x,*x2=0 


Figure  2-6.  Types  of  Deformation,  Defining  Piezoelectric 
Constants  d  and  e  (22:378). 

The  other  related  stress  constant,  e,  defines  another 
relationship  between  the  electric  field  and  the  stress  and 
strain.  Figure  2-6  gives  the  definition  of  e  graphically,  and 
the  equations  forms  are  shown  below.  That  is, 


( 2- 2 1  a) 


e_  mechanical  stress  developed  _/  dX  \ 
applied  electrical  field  ~\dElx 

c  _  el  ectrlcal  charge  densl  ty  developed  -idD  \  . 

applied  strain  \3x/e  '  ' 

Pyroelectric  Constant.  Some  materials  that  are 
piezoelectric  also  exhibit  pyroelectric  behavior. 
Pyroelectric  materials  absorb  thermal  energy  and  expand,  and 
the  reverse.  The  dimensional  changes  result  in  charge  being 
accumulated  (5:19).  The  current  that  is  developed  is 
proportional  to  the  temperature  change  rate.  The  pyroelectric 
constant  p  relates  the  charge  developed  (Q)  to  the  temperature 
change  (AT)  and  the  material's  cross-sectional  area  (A).  This 
relationship  is  defined  in  equation  form  as: 

Q~p*AT*A.  (2-22) 


The  voltage  (V)  due  to  the  pyroelectric  effect  is 
proportional  to  the  charge  developed  (Q)  divided  by  the 
material's  capacitance  (C)  (5:19).  That  is, 

v-  Q-pt(AT)A.pt(Ar)  (2_23) 

c  tA  e 

The  variables  represent  the  same  quantities  as  defined  in  the 
previous  examples. 

Piezoelectric  Materials 

There  are  two  main  groups  in  which  piezoelectric  materials 
can  be  classified:  natural  (quartz,  Rochelle  salt)  and 
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synthetic.  The  synthetic  group  can  be  further  divided  into 
crystals  (lithium  sulphur,  ammonium  dihydrogen  phosphate)  and 
polarized  f erroe 1 ectr i cs .  The  polarized  f erroe 1 ectr i cs  are 
typically  a  ceramic,  such  as  barium  titanate,  or  an  organic 
polymer  like  po 1 y v i ny 1 i dene  fluoride  (PVDF). 

Roche  lie  Salt.  Rochelle  salt,  (sodium  potassium 

tartrate),  which  has  a  chemical  formula  of  NaKC4H4Ot  •  4H,0 ,  has 
an  extremely  high  coupling  factor,  k,  along  its  f er roe  1 ectr i c 
axis,  e,  (17:17).  The  piezoelectric  efficiency  of  Rochelle 
salt  is  only  useful  between  its  lower  and  upper  Curie  points 
(255°  to  297*  K)  (17:17).  Outside  these  temperatures,  the 
piezoelectric  effect  is  greatly  reduced.  However,  as  pressure 
is  applied  to  Rochelle  salt,  its  upper  and  lower  Curie  points 
increase  linearly  (19:28). 

Rochelle  salt  is  environmentally  sensitive,  since  it 

contains  a  significant  amount  of  water  (18:115).  Mason,  in 

his  book  on  ultrasonics,  discusses  Rochelle  salt’s 

environmental  limitations.  He  says  that, 

...if  the  humidity  of  the  surrounding  atmosphere  is 
below  35  per  cent  at  25’  C,  the  water  vapor  pressure 
of  the  crystal  is  greater  than  the  vapor  pressure 
of  water  surrounding  atmosphere  and  the  crystal 
will  lose  water  and  dehydr ate ...  Above  85  per  cent 
humidity,  the  crystal  will  absorb  water  from  the 
atmosphere  on  its  surface  and  will  slowly  dissolve 
if  kept  is  such  an  atmosphere.  (18:117) 

Thus,  the  uses  of  Rochelle  salt  are  very  limited 

environmentally.  It  finds  its  primary  use  as  a  low  frequency 

transducer  (18:114). 
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Quartz .  Quartz  has  a  low  efficiency  of  mechanical  to 
electrical  conversion  ( e 1 ectr omechanc  i  a  1  coupling  factor,  k! )  . 
It  also  produces  only  a  small  charge  for  a  given  applied  load 
(9:1758).  Despite  a  low  coupling  factor,  quartz  has  found  use 
as  an  excellent  oscillator  due  to  its  very  stable  resonant 
frequency  (9:1758;  18:78). 

Cer ami cs .  A  ceramic  consists  of  a  conglomeration  of 
small  crystals  with  randomly  oriented  axes,  and  no  inherent 
piezoelectric  properties.  Fer roe  1 ectr  i  c  ceramics  must  be 
artificially  given  the  piezoelectric  effect  through  'poling' 
(  1  2:267  )  . 

This  process  can  be  described  as: 

...  the  polarity  needed  to  impart  p i ezoe 1 ectr i c 
properties  can  be  given  to  an  originally  isotropic 
polycrystalline  ceramic,  more  or  less  permanently, 
by  temporary  application  of  a  strong  electric 
field.  This  process,  called  "poling”,  is  analogous 
to  the  magnetizing  of  a  permanent  magnet.  (15:1) 

Typically  'poling'  is  done  at  high  voltages  (about  10  KV) 

while  the  material  is  heated  above  its  Curie  point,  then 

allowed  to  cool.  This  procedure  permanently  orients  the 

dipoles  in  the  same  direction,  making  it  piezoelectric 

(12:267).  Two  of  the  more  common  piezoelectric  ceramics  are 

barium  titanate  (Figure  2-1)  and  lead  zirconate  titanate 

(PZT)  . 

Organic  polymers.  Organic  polymers  as  piezoelectric 
materials  resulted  from  research  by  Kawai  of  Japan  (21:975). 
He  found  that  "the  polarized  homo  polymer  of  vinyl idene 
fluoride  (PVDF)  developed  far  greater  piezoelectric  activity 
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compared  to  any  other  synthetic  or  natural  polymer"  (5:1). 
PVDF  and  other  polymers  must  first  be  poled  before  they 
exhibit  p i ezoe 1 ectr i c  behavior. 


Pol i ng .  Polymers  can  be  poled  using  two  techniques: 
1)  application  of  a  large  DC  electric  field  or  2)  mechanical 
drawing.  The  electric  field  method  is  the  same  procedure  that 
is  used  with  ceramics  (see  Figure  2-7).  The  mechanical 
drawing  reorients  the  polymer  molecules  so  that  its  dipoles 
are  uniaxial  oriented,  but  not  necessarily  all  pointing  the 
same  direction  (22:326). 


Piezoelectric  Film  Dipole  Alignment  During  Poling 


Figure  2-7.  Piezoelectric  Film  Dipole  Alignment  During  the 
Poling  Process  (5:5). 


PVDF  Fabrication.  Of  the  many  available  polymers 
(such  as  vinylidene  fluoride  (VDF),  vinylidene  cyanide  (VDCN), 
and  po 1 y v i ny 1 i dene  fluoride  (PVDF)),  the  most  popular  one  is 
PVDF.  PVDF  is  a  sem  i -cr  y  sta  1  1  i  ne  polymer  which  has  long 
chains  of  repeating  units  of  CHj-CF,  (5:5).  In  this  (monomer) 
form,  more  than  90$  of  the  CH,-CF,  molecules  are  in  a  head-to- 
tail  form,  which  accounts  for  its  high  dipole  moment  (5:5). 
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Crystalline  forms  of  polyvinylidlne  fluoride 


i 


Figure  2-8.  Crystal  1 ine  Form  of  PVDF  Before 
and  After  Poling  (23-3). 


There  are  three  crystalline  phases  of  PVDF:  1)  alpha,  2) 
beta,  and  3)  gamma.  The  alpha-phase  results  after  the  polymer 
has  cooled  from  the  melt  (see  Figure  2-8);  this  phase  is 
non-polar  (5:5).  To  make  PVDF  piezoelectric,  it  first  must  be 
stretched.  This  stretching,  at  temperatures  below  the  melting 
point,  gives  rise  to  the  beta  form  which  has  carbon  chains  in 
parallel  strips  and  planes  with  pairs  of  hydrogen  and  fluorine 
atoms  randomly  distributed  about  the  carbon  chain  axis  (23:3). 
Before  poling  (which  gives  the  gamma-phase),  both  surfaces  of 
the  beta-phase  film  must  be  metallized  (see  Figure  2-9)  to 
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Active  axes  in  a 
piezoelectric  film 


3 


Figure  2-9.  Metallized  Film  Surfaces  for 
Poling  and  Film  Axes. 

create  electrodes  to  aide  in  poling  (23:3).  In  poling  (see 
Figure  2-7),  the  hydrogen  atoms  rotate  to  a  permanent  position 
toward  the  high-potential  side  of  the  carbon  atoms;  and  the 
fluorine  atoms  rotate  to  the  low  potential  side  of  the  carbon 
(see  Figure  2-8)  (23:3). 

PVDF  Pressure  Transducers 

A  piezoelectric  film  pressure  transducer  has  been 
developed  and  used  in  dynamic  situations  (3).  It  is 
inherently  smaller  and  cheaper  than  traditional  pressure 
transducers  which  use  piezoelectric  crystals  or  strain  gage 
bridges.  In  general,  it  possesses  the  following  advantages 
relative  to  conventional  piezoelectric  materials  (5:9): 

1.  It  has  a  broad  frequency  response  (DC  to  10  Mhz). 

2.  The  film  has  a  vast  dynamic  response  ( 10  1  to  1 0*  psi) 
so  it  could  detect  almost  any  pressure. 
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3.  The  high  electrical  impedance  of  the  film  (greater  than 
4X10I!  ohms)  closely  matches  the  impedance  of  CMOS 
devices,  which  reduces  the  electrical  current  leakage. 

4.  The  film  is  thin  and  flexible  so  that  it  can  be  cut, 
formed,  and  bonded  to  complex  or  simple  surface  shapes. 

5.  The  high  mechanical  strength  of  the  film  allows  it  to 
be  used  in  extreme  environments,  such  as  shock 
transition  zones. 

6.  Because  of  its  high  elastic  compliance,  it  accurately 
responds  to  pressure. 

7.  The  costs  to  fabricate  and  produce  a  film  transducer  are 
relatively  low. 

These  char acter i st i cs  of  piezoelectric  film  make  it  a  viable 
alternative  for  use  as  a  pressure  transducer. 

In  spite  of  its  many  desirable  char acter i st i cs , 
piezoelectric  film  does  possess  some  inherent  disadvantages  as 
a  pressure  transducer.  Since  piezoelectric,  film  is 
capacitive,  its  response  to  a  pressure  decays,  which  make  low 
frequency  measurements  difficult.  Piezoelectric  film  is 
pyroelectric  so  it  is  sensitive  to  temperature  changes.  The 
low  Curie  point  of  the  film  limits  its  use  to  below  100°  C. 
Any  surface  stresses  induced  by  the  surface  the  film  was 
mounted  to,  would  be  transmitted  to  the  film  as  well  as 
external  pressure  stresses.  Exposed  electrodes  are  sensitive 
to  electromagnetic  radiation.  This  sensitivity  would  add 
error  to  a  pressure  measurement.  Despite  these  limitations, 
most  of  the  film’s  ill  effects  can  be  minimized  or  eliminated. 

Summary 

The  piezoelectric  effect  has  long  been  known,  but  not 
under  stood  unt  i  1  recent  times.  This  property  occurs  naturally 
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in  some  crystals,  and  it  can  be  artificially  induced  in 
materials  such  as  ceramics  and  polymers  through  'poling'. 
P i ezoe 1 ectr i c  materials  exhibit  many  properties,  in  the  case 
of  a  pressure  transducer,  the  ability  to  convert  a  force  to  an 
electrical  signal  makes  piezoelectric  transducers  possible. 
Of  the  many  piezoelectric  materials,  PVDF  film  has  the 
potential  to  be  very  useful  for  pressure  measurements.  It 
would  be  small  and  relatively  inexpensive  compared  to 
conventional  p i ezoe 1 ectr i c  transducers.  Though  PVDF  does  have 
limitations  as  a  pressure  sensor,  most  of  these  can  be  reduced 
or  eliminated  to  make  it  a  promising  new  type  of  transducer. 
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III.  Experimental  Design.  Development,  and  Evaluation 

Procedur es 


To  use  piezoelectric  film  as  a  low-frequency  and  dynamic 
pressure  transducer,  several  techniques  must  be  developed 
( samp  1 i ng- and-reset  and  time-history  reconstruction).  This 
task  involves  designing  and  developing  an  analog  circuit  to 
discharge  the  accumulated  charge  of  the  piezoelectric  film. 
Also,  a  technique  to  identify  reset  pulses  and  to  accurately 
reconstruct  the  data  must  be  developed. 

Reset  Circuit  Design 

There  are  several  aspects  that  must  be  considered  in 
order  to  properly  design  a  reset  circuit.  The  equivalent 
electrical  circuit  of  the  film  needs  to  be  established  to 
ensure  an  accurate  interface.  Since  piezoelectric  film 
manifests  electrical  capacitance  and  its  charge  decays,  this 
property  needs  to  be  minimized.  The  reset  circuit  must  have 
a  response  sufficiently  broad  to  respond  to  the  desired 
frequency  range.  Electrical  components  that  are  interfaced 
must  have  matched  impedances  to  minimize  charge  leakage  and 
sources  of  electrical  noise. 

Piezoelectric  Film  Equivalent  Circuit.  Piezoelectric 
film  can  be  represented  as  a  charge  or  current  generator  (as 
shown  in  Figure  3-1)  with  film  capacitance,  Cf ,  and  film 
resistance,  Rf  (5:27).  The  p i ezoe 1 ectr  i  c  f i lm  capac i tance  ,  C, , 
is  proportional  to  the  electrode  area  and  inversely 
proportional  to  the  f  i  lm '  s  th  ickness  (  5  : 27  ) .  The  charge  generated 
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Figure  3-1.  Equivalent  Electrical 
Circuit  for  the  Piezoelectric  Film  as 
a  Charge  Generator  (5:27). 

by  the  film  is  proportional  to  the  applied  force  or  induced 
displacement,  and  it  can  be  expressed  as: 

Q~dt*F-et*xt  (3-1) 

where  dt  and  et  are  piezoelectric  constants,  F  is  the  force, 
and  xt  is  the  deflection  along  the  film’s  thickness  (5:27; 
12:269).  In  order  to  consider  a  more  useful  form,  the  charge 
generator  can  be  converted  to  a  current  generator  as 
illustrated  in  equation  3-2  and  depicted  in  Figures  3-1  and 


(3-2) 


3-2.  The  current  generated  by  the  film  is  the  sum  of  the 
currents  flowing  in  Cf  and  Rf ;  that  is, 


+  i 


(3-3) 


Since  the  output  voltage  (V,)  is  the  ratio  of  the  charge  (Q) 
to  the  film  capacitance  (Cr),  it  can  be  written  as: 


^(/V dt)/C'=(f(ic-l,,)dt)/cr  (3-4) 
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Figure  3-2.  Equivalent  Electrical  Circuit  for  the 

Piezoelectric  Film  as  a  Current  Generator  (12:270). 

The  above  equation  can  be  solved  for  (if  -  ilf)  by  making 
substitutions  for  ir  and  i,f.  That  is, 


if 


(3-5) 


Rearranging  terms  and  solving  for  the  ratio  of  the  output 
voltage  to  the  displacement  gives: 


KxD 
xD  +  1 


(3-6) 


where  tf  A  time  constant  =  Rf*Cf  (sec) 

K  A  sensitivity  =  et/Cf  (V/cm) 

D  A  derivative  =  d  /dt  (12:270). 
Equation  3-6  can  be  rewritten  to  establish  the  time  response 
of  the  piezoelectric  film  (12:270).  That  is, 


(xD  +  l)V0-(xKD)Xt  .  (3-7) 

If  an  external  force  is  suddenly  applied  and  allowed  to  remain 
on  a  piece  of  PVDF  film,  it  would  cause  the  displacement  "A" 
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Figure  3-3.  Charge  Decay 
of  a  Piezoelectric  Device 
(  12:271  )  . 

and  the  voltage  increase  shown  in  Figure  3-3.  Since  "A"  is 
constant  for  the  time  duration  T,  that  is,  0  <  t  <  T,  the 
right-hand  side  of  equation  3-7  is  zero.  That  is, 

(tD  +  1)Vo  =  0  .  (3-8) 

At  t=0,  when  the  force  was  applied,  the  charge  instantaneously 
increased  to  et*A,  and  it  gives  the  following  initial 
condition  (12:271): 

V0-^-  at  t-0 *  .  (3-9) 

ct 

Solving  equation  3-8  with  the  initial  condition  gives  the 
following  result: 

V0-(etA/Ct)  exp(-t/tf)  for  0<t<T.  (3-10) 

At  t  =  T,  the  force  is  removed  and  returns  to  its  original 
value.  This  displacement  behavior  causes  the  same  magnitude 
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of  response  voltage  to  be  observed,  but  it  is  of  opposite 
polarity.  Thus,  V.  decreases  by  et*A/C,  at  t  =  T  (12:271).  As 
a  result  at  t  =  T' ,  V,  is  expressed  as 

vo--^  (exp|--£j  -lj  exp|^tc~-T)|  for  t>T  .  (3-11) 

This  is  the  initial  condition  for  t  >  T.  When  it  is 
substituted  into  equation  3-8  and  solved,  the  solution  becomes 

<3-,2) 

which  gives  the  solution  for  all  time  greater  than  T,  the 
force  release  time  (12:271).  As  the  time  constant  t  is 
increased,  the  response  to  x,  :  A  becomes  more  linear. 

When  the  film  is  connected  to  a  cable  and  some  type  of 
amplifier,  this  arrangement  changes  the  response  of  the  film. 
These  effects  must  be  accounted  for  in  the  reset  circuit 
design.  Figure  3-4  illustrates  the  typical  equivalent  circuit 
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Figure  3-4.  Equivalent  Electrical  Circuit  for  the 

Piezoelectric  Film,  Cable,  and  Amplifier  (12:270). 
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for  a  piezoelectric  film  measuring  circuit,  cable,  and 


amplifier  circuit.  These  additional  electrical  impedance 
effects  can  be  easily  combined  with  the  film’s  capacitance  and 
resistance,  as  shown  in  equations  3-13  and  3-14.  That  is, 


R _ +  Rf 


(3-13) 


C  -  Cf  +  CcaJ>la  +  Cam> 


(3-14) 


The  cable  resistance  is  neglected  since  it  is  small  compared 
to  the  film  or  amplifier's  resistance.  Overall,  these  new 
values  (R  and  C)  need  to  be  substituted  into  the  equations  for 
the  film,  to  account  for  their  cumulative  effect. 

Circuit  Design.  There  are  two  fundamental  types  of 
analog  circuits  that  could  be  used  to  measure  the 
piezoelectric  film's  response.  A  voltage  follower  circuit,  as 
shown  in  Figure  3-5,  and  a  charge  amplifier  or  integrator 
represented  by  Figure  3-6.  Each  circuit  has  distinct 
advantages  and  disadvantages. 


Reset  swi  tch 


Sampled  output 


PVDF 

Film 


Sampling  command 


Figure  3-5.  Voltage  Follower  Circuit  (14:27). 
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The  voltage  follower  circuit  (Figure  3-5)  is  very  simple 
and  it  can  easily  be  modified  to  achieve  different  levels  of 
gain.  It  is  very  sensitive  to  capacitive  inputs,  such  as  from 
the  film  or  cable  (14:105).  A  very  high  input  impedance 
amplifier  is  needed  to  reduce  charge  leakage,  which  in  turn, 
increases  noise  sensitivity.  The  rate  of  force  application 
speed  has  an  affect  on  the  amplifier’s  output,  since  typical 
high  impedance  amplifiers  have  a  slow  slew  rate  (14:106). 


The  charge  amplifier  circuit  (Figure  3-6)  is  not 
influenced  by  the  time  constants  of  the  piezoelectric  film  or 
the  cable  due  to  their  capacitive  nature.  Since  a  charge 
amplifier  is  current  driven  with  a  zero  input  impedance,  no 
voltage  is  generated  across  the  film  (5:30).  Because  the 
charge  amplifier  is  immune  to  capacitance  influences,  looser 
tolerances  can  be  used  on  film  size  and  cable  length  (5:30, 
12:718).  Nor  is  the  amplifier  affected  by  variations  of  the 
film's  dielectric  constant  due  to  temperature  change  (14:105). 


In  summary,  the  charge  amplifier  circuit  possesses  an 
abundance  of  desirable  performance  properties  for  capacitive 
type  transducers.  Thus,  a  charge  amplifier  is  the  circuit  of 
cho ice. 

Sampling  and  Reset  Operation 

The  conventional  charge  amplifier  would  instantaneously 
indicate  the  force  which  is  applied  to  its  associated  sensor. 
This  behavior  can  cause  problems  if  the  input  force  is  either 
a  low  frequency  excitation  or  it  has  a  large  magnitude. 

The  charge  amplifier  cannot  accurately  read  low  frequencies 
from  a  piezoelectric  sensor  since  a  piezoelectric  input  decays 
with  respect  to  time.  The  charge  decay  could  be  larger  than 
the  input  force  response.  If  the  input  force  is  highly 
dynamic,  the  input  could  saturate  the  amplifier's  input,  which 
would  result  in  a  saturated  output  (7). 

The  sampling  and  reset  operational  mode  overcomes  these 
limitations,  and  it  is  illustrated  in  Figure  3-7.  The  circuit 
implementation  is  shown  with  a  charge  amplifier  in  Figure  3-8. 
When  the  reset  switch  is  open,  the  charge  amplifier  follows 
the  input  from  the  piezoelectric  device,  and  the  output  is 
sampled  as  shown  in  Figure  3-7d.  The  input  variations  during 
the  sampling  period  appear  as  an  amplitude  change  of  the 
sampled  output  (7:1).  After  sampling,  the  switch  is  closed  to 
discharge  the  charge  accumulation  on  the  amplifier's 
capacitor.  Again,  after  the  switch  opens,  the  amplifier 
follows  the  change  of  the  input  signal  (7:2).  This  process 
practical ly  el iminates  amp  1  if ier  saturation,  which  means  that 
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Figure  3-7.  Sampling  and  Reset 
Operation  (14:25). 


the  amplifier  can  operate  with  a  broad  and  dynamic  response. 
Low  frequency  measurements  become  possible  since  charge  decay 
is  minimized. 

Reset  Implementation.  Figure  3-8  depicts  the  situation 
where  the  piezoelectric  film  is  modelled  as  a  capacitor  Cf, 
which  has  an  output  proportional  to  an  applied  force  (14:24). 
The  reset  switch  could  be  mechanical  or  electrical,  but  in 
this  development,  it  is  electrical.  It  is  character i zed  as 
having  an  R„  of  about  400  ohms,  and  an  Rt(f  of  IE  15  ohms,  or 
virtually  an  open  circuit. 
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Charge  Amplifier  Reset  and  Sampling  Scheme 
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Figure  3-8. 

(14:26). 
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Reset  State .  When  the  reset  switch  is  closed,  it  is 


modelled  as  R91  in  Figure  3-8b.  The  closed  switch  causes  the 
film's  charge  to  be  transferred  to  the  feedback  capacitor  of 
the  charge  amplifier,  and  it  is  discharged  through  Rori  in  a 
short  period  of  time  called  TPtttl  (14:24)  Since  the  negative 
input  terminal  of  an  operational  amplifier  is  at  virtual 
ground  level,  the  charge  is  effectively  grounded  through  the 
reset  technique  (see  Figure  3-7d)  (13:178).  During  the  reset 
period,  the  charge  cannot  be  measured,  even  though  the  force 
may  still  be  applied;  thus,  the  reset  command  pulse  width  Trmt 
(see  Figure  3-9)  should  be  minimized  (14:28).  The  minimum 
reset  time  is  a  function  of  the  time  constant  of  the  reset 
circuit  tr„,t  (14:28).  That  is, 

riM«t  >  xt0M0t  =  Roa  *  Cfto  •  (3-15) 

Typically,  the  reset  pulse  width  should  be  established  for  3 
to  5  time  constants  (tf„lt)  so  that  most  of  the  charge  can  be 
dissipated.  If  T„„t  =  95%  of  the  charge  will  be 

dissipated,  and  if  TrMlt  =  5*trMtt,  it  will  dissipate  99.3%  of 
the  charge  (15:9). 

Samp  1 i ng  State .  As  the  reset  pulse  ends  (goes  low) 
and  the  switch  in  Figure  3-8c  opens,  the  amplifier  follows  the 
input  signal,  as  illustrated  in  Figure  3-7d  (14:29).  The 

signal  flows  predominantly  through  the  capacitor,  since  R,ff  is 
like  an  open  circuit.  Now,  the  amplifier's  output  signal  can 
be  sampled  by  the  A/D  converter.  The  sampling  period  lasts 
for  the  time  duration  (t,  -  t,.,).  Next,  the  reset  pulse  goes 
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high,  and  the  reset  process  repeats.  To  get  the  actual  time 
history,  the  sampled  outputs  must  be  integrated  to  recreate 
the  equivalent  output  history. 

Reset  Circuit  Characteristics 

Offset  Calibration.  The  reset  scheme  eliminates  the  need 
to  do  a  zero  offset  calibration  for  the  measurement.  Each 
reset  cycle  clears  the  accumulated  charge  due  to  any  initial 
conditions  (14:31).  However,  the  operational  amplifier  must 
be  adjusted  for  zero  offset,  but  this  situation  is  independent 
of  the  PVDF  film. 

Sensitivity  and  Dynamic  Range.  The  conventional  charge 
amplifier  is  limited  by  sensitivity  and  dynamic  range  (14:31). 
However,  the  product  of  the  two  is  a  constant.  That  is, 

Sens i  tivi  tyX  Dynamic  Range  -  Constant  .  (3-16) 

The  constant  is  determined  by  the  saturation  limits  of  the 
charge  amplifier;  it  is  typically  ±15V.  Thus,  the  sensitivity 
can  be  increased  by  a  decrease  in  the  dynamic  range,  and  vice 
versa  (14:31). 

The  reset  scheme  minimizes  this  interdependence,  since  a 
higher  reset  frequency  allows  less  of  an  input  amplitude 
increase  before  it  is  reset.  Thus,  the  interdependency  is 
reduced . 

The  limit  of  the  sensitivity  associated  with  the  reset 
scheme  is  bounded  by  the  saturation  voltage  of  the  charge 
amplifier  (Vut).  This  relationship  is  expressed  below  by  the 
saturation  voltage  limit.  That  is, 
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I*pl  (f)4r=lspl  (f)(7)  <  I' "’-I 


(3-17) 


or 

\Sp^\  =  \Vaat\/{^{^  (3-18) 

where  S,  A  (dl*X3*t)/e  =  film  sensitivity  (V/lb,) 

t  A  film  thickness  (pm) 

£  A  film  dielectric  constant  (F/m) 
dF/dt  A  force  application  speed  (lbf/sec) 
f  A  reset  command  rate  (Hz) 

Vllt  A  charge  amplifier  saturation  voltage  (V). 
As  a  result,  a  higher  reset  rate  or  slower  force  application 
speed  increases  the  sensitivity  (14:32,33). 

Low  Frequency  Response.  As  mentioned  previously, 
piezoelectric  film  cannot  measure  static  forces  since  its 
charge  decays  with  time.  With  the  reset  process,  the 
developed  charge  is  sampled  and  discharged  before  the  charge 
decay  process  manifests  a  significant  effect.  Thus,  the  reset 
rate  must  be  sufficiently  short  for  the  amplifier  to  operate 
in  the  "linear"  portion  of  the  charge  curve.  For  practical 
purposes,  this  is  limited  to  about  one-fifth  of  the  circuit 
time  constant  (15:14).  Thus,  a  larger  circuit  time  constant 
permits  a  slower  reset  rate. 

To  better  understand  the  effect  of  the  reset  frequency  on 
minimizing  the  charge  decay,  a  linear  input  force  application 
process  can  be  examined.  The  applied  force  causes  a  charge  to 
be  developed  on  the  film  that  is  proportional  to  the  force 
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application  rate, 
f o 1  1 ows  : 


This  situation  can  be  represented  as 


c-  dQcJt  -  charge  Development  Rate  (  —  °—  \  (3-19) 

dt  \  sec  / 

where  c  =  Cf  *  S,  *  dF/dt 

dF/dt  =  force  application  rate. 

Since  the  PVDF  film's  output  voltage  is  the  developed  charge 

divided  by  the  film  capacitance,  the  film's  residual  charge 

can  be  found  by  multiplying  equation  3-10  by  Cf  to  get: 

0<t)  =0oexp(-^j  (3-20) 

where  Q,  =  et*A  =  an  initial  charge. 

But  while  the  charge  is  accumulating,  it  is  also  decaying  by 
an  amount,  Qt .  Multiplying  equation  3-11  by  C,  gives  an 
expression  for  Q, . 

Qa  ( t)  =0o|l-exp|-^jj  .  (3-21) 

The  difference  between  the  charging  and  discharging  rates  is 
shown  below.  That  is, 

dQ(t)  .  dQc(E>  _  d0d^  (3-22) 

dt  dt  dt 

Solving  equation  3-22  yields  the  observable  charge  at  an 
instant  of  time  t.  That  is,  when  t  =  T,  the  charge  is: 
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Q(T)  =x  c|l- exp . 


(3-23) 


If  t  were  infinite,  the  output  would  be  linear;  and  the 
observed  charge  would  be  equal  to  the  developed  charge 
(14:36).  Since  t  is  indeed  finite,  at  best,  the  observed 
charge  will  be  t*c. 

With  the  reset  scheme,  the  time  T  becomes  AT,  the  time 
between  reset  cycles.  The  accumulated  charge  for  each 
sampling  period  becomes: 

0(AD  =xc|l-exp(-lAlJJ  .  (3-24) 

The  complete  charge  time  history  at  time  T  is  given  by  the 
summation  of  the  sampled  outputs  (14:38).  That  is, 

0(T)  =  J2  (x  c)  (1  -exp(ATA)) 

=(^)tC(1-eXP(~Ar/t))  '  (3’25) 

-  x  c  T  f  (1  -  exp(-  1/t  f )) 

If  equation  3-24  is  divided  by  cT,  the  result  is  the  ratio  of 
the  observed  charge  to  the  developed  charge.  By  letting  this 
ratio  define  an  error  function,  called  Q,  (14:38),  then 

O.-^r  -«*(•=$))  .  (3-26) 

As  the  product  x*f  is  increased,  the  charge  error  decreases. 
When  t*f  is  100,  the  charge  error  is  less  than  one  percent  (Q, 
=  0.9950),  and  each  decade  improvement  of  t*f  decreases  the 
charge  error  by  a  similar  amount  (14:40).  Table  3-1  shows  the 
charge  error  for  several  i*f  products. 
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Table  3-1.  Charge  Error  for  Various  t*f  Products  (14:40). 


Q.  (%) 

0.01 

1  .  00 

0  .  1 

10.00 

1 

63.21 

1  3 

85.04 

5 

90 .63 

10 

95.16 

100 

99.50  | 

1000 

99.95 

10000 

99.995 

This  behavior  shows  that  Q#  is  minimized  by  a  large  t*f 
product.  In  this  case,  a  shorter  reset  rate  is  required  for 
a  larger  t ,  or  as  r  decreases,  f  must  increase  to  maintain  the 
same  level  of  accuracy  (14:39). 

Reset  Rate  Limits.  The  maximum  reset  rate  is  determined 
by  the  minimum  detectable  signal  of  the  charge  amplifier  or 
the  minimum  resolution  of  the  A/D  device.  For  a  particular 
A/D  device,  the  minimum  level  is  defined  by  the  resolution  of 
the  Least  Significant  Bit  ( LSB )  .  From  equation  3-16,  the 
minimum  detectable  voltage  is 

l^«J  =  |Spl  \dF/dt\(l/f)  (3-16) 
or 

.  (3-27) 

Hence,  the  upper  reset  rate  can  be  increased  by:  1)  increasing 
the  sensitivity;  2)  lowering  the  minimum  detectable  signal; 
or  3)  increasing  the  force  application  rate  (14:48). 
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The  charge  decay  time  constant  (ttifCljt)  and  charge 
development  rate  can  also  determine  the  upper  reset  frequency. 
From  equation  3-23 


Q(AT)  =tc(l-exp(^i))=|0]lin|  .  (3-23) 

When  Q(AT)  =  |  Q.j, |  =  |v,j,|*c>  the  maximum  reset  rate  can  be 

obtained  from  equation  3-28  (14:48).  That  is, 

ln[l  '  (3‘28) 

The  lower  reset  frequency  is  likewise  determined  by  the 
charge  amplifier  limits.  In  this  case,  the  voltage  saturation 
level  (VIlt)  sets  this  lower  bound.  Equation  3-16  can  be  used 
to  determine  the  saturation  voltage.  That  is, 


\vMt\>  \sp\\dF/dt\a/f)  . 


(3-16) 


Solving  for  f  yields  the  following: 


(3-29) 


The  circuit  time  constant  and  charge  development  rate  can  also 
determine  the  lower  reset  rate.  Similar  to  the  development 
for  the  lower  rate  can  be  determined  by  equation  3-30 
(14:49).  That  is, 


‘  rmt.!. ' 


_ 1 _ 

l^l  ~(Ccixcult\V*at\  /  xciTcultc)\ 


(3-30) 
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Circuit  Implementation 


Two  of  the  primary  purposes  of  the  charge  amplifier  with 
reset  are  to  minimize  charge  leakage,  and  then  to  discharge 
the  accumu lated  film  charge  during  the  reset  period.  Ideally, 
the  input  impedance  of  the  operational  amplifier  would  be 
infinite  and  would  draw  no  current  or  have  a  zero  input  bias 
current.  Similarly,  the  reset  switch  would  have  a  zero  R,„ 
resistance  and  infinite  Rsff  resistance.  Since  neither  of 
these  conditions  are  possible,  a  circuit  that  comes  close  to 
satisfying  these  conditions  will  suffice.  If  the  simple 
resetting  charge  amplifier  in  Figure  3-6  were  to  be  used  with 
a  FET  reset  switch,  the  dra i n-to-source  current  leakage  would 
cause  a  small  current  to  flow  into  the  summing  junction  when 
the  switch  is  off.  This  situation  would  dissipate  the  charge 
on  the  film  and  increase  the  measurement  error  (13:223).  The 
charge  amplifier  circuit,  with  reset  configuration  shown  in 
Figure  3-11,  minimizes  the  current  leakage  and  has  a  low  R„ 
resistance.  Even  with  an  ultra-low  bias  current  operational 
amplifier  and  a  low  leakage  feedback  or  integrating  capacitor, 
the  summing  junction  current  leakage  can  be  the  dominant  error 
source  (13:224).  The  reset  part  of  the  circuit  uses  two 
n-channel,  enhancement-mode  MOSFETs  (Q,  and  Qt )  as  the  reset 
switch.  Both  MOSFETs  are  switched  simultaneously.  Q,  is 
switched  with  gate  voltages  of  zero  and  plus  ten  volts.  This 
situation  eliminates  gate  and  source-to-drain  current  leakage 
during  the  sampling  period  when  the  switches  are  off 
(  zero  gate  vo  1  tage)  and  the  reset  command  pulse  is  low  (13:224). 
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Figure  3-10.  Charge  Amplifier 
with  a  Low-Leakage  Reset  Switch 
( 73:224) . 


When  the  reset  pulse  command  is  high,  the  MOSFET  switches 
turn-on,  and  the  capacitor  is  discharged  with  the  effective 
resistance  R„  ,rf .  The  R„  ,ff  resistance  is  the  sum  of  R„  for  Q, 
and  Qt  =  200Q  +  200Q  =  400Q).  When  the  reset  command 
goes  low  again  and  the  MOSFETs  are  off,  the  small  leakage 
current  through  Q,  passes  through  R,  to  ground  with  a 
negl igible  voltage  drop  (13:224).  The  summing  junction  has  no 
current  leakage  since  the  source,  drain,  and  substrate  of  Q,, 
are  all  at  the  same  voltage  level.  Thus,  current  leakage 
during  the  off  cycle  is  minimized  through  the  reset  switch. 


3-20 


Component  Selection.  With  the  proper  selecting 

components,  the  current  leakage  can  be  reduced.  Typical 
operational  amplifiers  have  an  input  impedance  greater  than 
2  MQ  and  an  associated  input  bias  current  which  is  typically 
several  hundred  nanoamps  in  magnitude.  The  large  impedance 
mismatch  between  the  film  (Rf  >  4.2E12  Q)  and  the  operational 
amplifier,  along  with  the  relatively  large  input  current, 
would  dissipate  charge  from  the  film.  An  ultra-low  bias 
current  JFET  operational  amplifier,  such  as  the  OPA128  which 
is  typically  called  an  electrometer  operational  amplifier,  has 
an  input  impedance  of  1E15  ohms  and  an  input  bias  current  as 
low  as  60E-15  amps  (60  femtoamps)  (25:2-64).  This  type  of 
operational  amplifier  has  an  input  impedance  greater  than  the 
film's  impedance,  and  this  greatly  reduces  current  leakage  as 
compared  to  the  "typical"  operational  amplifier. 

The  integrating  or  feedback  capacitor  can  contribute  to 
current  leakage.  Standard  film  capacitors  have  an  insulation 
resistance  (IR)  (resistance  to  a  DC  current  flow  under  steady 
state  conditions)  on  the  order  of  50,000  to  100,000  MQ-pF.  A 
high  quality  teflon  or  polystyrene  capacitor  has  an  IR  greater 
than  1,000,000  MQ-pF.  Assuming  a  0.02  pF  feedback  capacitor 
is  used,  the  capacitor's  leakage  resistance  is  on  the  order  of 
5E7  MQ  (5E13  ohms)  (14:224). 

If  components  such  as  these  are  used,  the  circuit  leakage 
currents  (not  including  the  reset  switch)  at  the  summing 
junction,  under  a  worse  case  of  10V  full  scale,  are  less  than 
1  picoamp  (14:224).  A  common  switching  MOSFET,  such  as  the 
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3N170,  has  a  maximum  specified  leakage  current  of  10  nanoamps. 
As  a  result,  the  MOSFET  switch  contributes  up  to  10,000  times 
as  much  leakage  current  as  does  the  capacitor  and  operational 
amplifier  combined  (14:224).  Thus,  the  overall  circuit  design 
is  very  important  in  order  to  minimize  the  MOSFET’s  source 
leakage  current.  Properly  selected  components  can  further 
reduce  this  current  leakage  to  a  tolerable  level. 

Design  Calculations 

Before  some  of  the  component  values  can  be  established, 
the  expected  film  properties  and  output  response  must  be 
known.  These  values  are  based  upon  the  anticipated  PVDF 
pressure  transducer's  operating  environment  and  equipment 
limitations.  Some  of  the  equipment  limitations  are  imposed 
locally,  and  not  by  the  commercial  availability  of  better 
equ i pment . 

The  equipment  limitations  include:  1)  32,000  Hz  sampling 
rate  limit  for  the  A/D  board;  2)  12-bit  resolution  of  the  A/D 
board;  3)  memory  limitation  of  4096  sampled  points;  and  4) 
maximum  reset  pulse  amplitude  of  ±10V.  These  limitations  will 
be  discussed  as  they  are  encountered  in  the  following 
d i scuss ion . 

The  design  was  based  on  the  following  operational  limits: 
1)  a  pressure  range  of  0-2  pounds  per  square  inch;  and  2)  a 
frequency  band  from  DC  to  16,000  Hz.  The  pressure  response  is 
typical  for  commercially  available  dynamic  pressure 
transducers.  It  is  not  imposed  by  the  film,  but  it  does 
affect  the  choice  of  secondary  gain  and  A/D  device  selection. 
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The  frequency  band  width  is  imposed  by  the  available  A/D 
board,  and  this  limit  could  easily  be  improved  with  a  faster 
sampling  board.  However,  with  these  constraints,  the  film 
properties  can  be  calculated. 

Film  Calculations.  The  film  size  was  selected  to  be  1 
square  centimeter.  This  feature  was  not  based  upon  any 
optimized  size,  but  rather,  it  was  influenced  by  the  ease  of 
handling,  fabrication,  and  a  standard  unit  of  length.  The 
voltage  output  for  the  film  under  a  pressure  field  is 
proportional  to  its  thickness  and  is  independent  of  its  area 
(5:21).  Thus,  smaller  film  areas  would  produce  the  same 
output  voltage,  but  a  smaller  amount  of  accumulated  charge, 
which  is  proportional  to  the  film's  area  (5:21). 

F  i  1m  Capac  i  tance  (Cf)  and  Sensitivity  (S().  The  film 
sensitivity  (S#)  indicates  the  output  voltage  per  unit  of 
applied  force.  Since  the  prospective  piezoelectric  film 
transducer  will  measure  pressure,  its  output  needs  to  have  the 
units  of  volts  per  square  centimeter.  Equation  3-31 
establishes  the  relationship  between  film  sensitivity,  film 
constants,  and  film  dimensions.  That  is, 

S~  *3  C  ( mv/psi )  .  (3-31) 

p  e 

Film  capacitance  is  also  an  important  parameter  involved  in 
many  calculations  and  it  is  found  from  the  following 
expression.  That  is, 
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C-eA/t  . 


(3-32) 


In  both  equations,  t  is  the  film's  thickness,  e  is  the  film's 
dielectric  constant  (106E-12  (Farads/meter ) ) ,  is  the 
applied  pressure,  A  is  the  electrode  area  (1  cm1),  and  dt  is 
the  piezoelectric  strain  constant  (-22E-12  (C/cml)/(N/ml) )  . 


The  values  for  Cf  and  S, 

Table  3- 2 . 

Table  3-2.  Film 

for  three  film  thickness  are  shown  in 

Capacitance  and  Sensitivity. 

c,  (df) 

28 

378.6 

40.10 

52 

203.8 

74.40 

110 

96.40 

157.4 

Reset  Rate  Limit  Calculation.  The  upper  and  lower  reset 
rates  were  calculated  from  the  equations  given  in  the  Reset 
Rate  Limits  section.  Before  this  calculation  can  be 
accomplished,  several  other  constants  must  be  found.  The 
minimum  detectable  voltage  is  determined  by  the  LSB  resolution 
of  the  A/D  board  which  is  5V/4095  =  1.22  mv .  The  force 
application  rate  is  limited  by  the  maximum  sampling  rate  of 
32,000  Hz.  Since  10  samples  per  reset  period  are  used  to 
accurately  reconstruct  the  time  history,  the  highest  force 
application  rate  would  correspond  to  2  psi  in  one-half  of  a 
samp  1 i ng  eye  1 e . 
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r  2  psi  *  ~?2'—00  —  ♦  0.5-12, 800  (3-33) 

dt  10  sec 

To  find  the  reset  limits  imposed  by  the  circuit's  time 

constant  and  charge  development  rate,  these  constants  (tcjrc|jl 

and  c)  need  to  be  calculated.  From  equations  3-13  and  3-14, 

tcireiit  can  be  calculated  based  upon  the  values  for  an  OPA128 

type  operational  amplifier,  Rf  and  Cf  for  each  film  thickness, 

and  a  Ccltll  of  50  pF .  For  28  pm  film  tcjrcijt  is  1,798  seconds, 

which  is  4.2  X  1 0 11  ohms  times  429  pF.  The  value  for  c  (charge 

development  rate)  is  determined  from  the  following  equation: 

.  (3-34) 

*  at 

Table  3-3  summarizes  the  results  for  c  and  tejr{lH  as  a  function 
of  film  th i ckness . 


Table  3-3.  Circuit  Time  Constant  and  Charge 
Development  Rate. 


Film  Th i ckness 

C, 

^circtit 

c 

Rf 

(pm) 

(PF) 

(s) 

(nC/s) 

(Q) 

28 

379 

"  "  . . .  ‘ 

1798 

52 

204 

1971 

194.3 

7.8E12 

110 

96.4 

- 1 

2398 

193.4 

16.5E12  | 

Now  that 

al  1 

the  required 

par ameter 

values 

have  been 

ca 1 cu 1 ated 

f... 

and  fijt  can  be 

ca 1 cu 1 ated . 

These 

values  are 

shown  in  Table  3-4  for  each  film  thickness. 
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Table  3-4.  Reset  Rate  Limits. 


The  minimum  reset  rate  in  each  case  is  determined  by  the 
operational  amplifier  saturation  voltage,  not  the  circuit 
constants.  The  theoretical  upper  limit  is  set  by  the  minimum 
detectable  voltage  of  the  A/D  board;  but  in  fact,  the  highest 
reset  rate  is  imposed  by  (1/10  of  32,000  Hz  or  3200  Hz)  the 
A/D  board.  Thus,  the  actual  reset  rate  can  range  from  51  Hz 
to  3200  Hz  depending  upon  the  film  thickness. 

An  initial  reset  rate  of  500  Hz  is  selected  so  that  the 
charge  error  will  be  less  than  0.12%.  This  condition  will 
allow  up  to  0.82  seconds  (4096  samples)  of  the  sampled  data  to 
be  recorded.  Although  this  time  duration  is  not  long,  it  is 
adequate  for  initial  testing.  In  an  actual  application,  more 
memory  would  be  required  to  support  a  longer  test  run  or  a 
higher  sampling  rate. 

Reset  Pulse  W_j  dth .  The  reset  pulse,  which  is  determined 
by  the  time  constant  of  the  reset  switch  and  feedback 
capacitor,  can  be  calculated  using  equation  3-14.  Using  a 
pair  of  3N170  MOSFETs  with  R„  =  200  ohms  each,  and  a  0.02  pF 
feedback  capacitor,  rrMll  becomes: 
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trM#t  -  0 . 02  *  2  *  200  Q  =  8  |isec  . 


(3-35) 


Assuming  Trtlll  =  5  *  rrillt  so  that  over  99%  of  the  charge  is 
dissipated,  Triltt  is  40  ps.  This  value  gives  a  starting  point 
for  circuit  testing. 

Time  History  Reconstruction 

The  piezoelectric  film's  output  is  reset  periodically  to 
eliminate  any  charge  decay  error.  Each  reset  period  is 
composed  of  ten  data  points.  To  get  the  entire  history,  the 
sampling  periods  have  to  be  integrated  and  the  error,  induced 
by  the  DC  offset  and  rise  rate,  must  be  removed. 

This  procedure  involves  locating  the  beginning  of  each 
reset  period  and  adding  that  point  to  the  last  point  of  the 
prior  reset  period  sample.  The  process  of  locating  each 
initial  portion  of  a  collected  sample  can  be  difficult  for 
low-level  signals.  There  are  many  ways  to  do  this,  but  two  of 
the  simpler  methods  are:  1)  to  synchronize  the  reset  pulse  and 
digital  sampling;  and  2)  to  make  the  reset  pulse  sufficiently 
wide  so  that  at  least  one  sample  data  point  is  taken  during 
the  reset  pulse  width,  giving  a  saturated  signal. 

The  second  method  was  chosen  because  it  was  simpler  and 
since  the  purpose  of  this  thesis  was  to  demonstrate  a  concept 
rather  than  accomplish  a  complete  sensor  development.  The 
pulse  width  was  selected  to  be  100  ps,  as  opposed  to  a  minimum 
of  approximately  40  ps  needed  to  discharge  the  feedback 
capac i tor  . 
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Once  the  sample  interval’s  beginning  and  end  are 
established,  the  first  point  of  the  new  reset  period  is 
erroneous  and  cannot  be  used  as  data.  Thus,  a  simple  straight 
line  interpolation  is  made  between  the  previous  end  point  and 
the  second  data  point  measured  with  respect  to  the  beginning 
of  the  data  record.  In  equation  form,  the  value  for  the 
"bridge"  point  can  be  found  as  follows: 


*n  = 


-  *a-l  +  (*a-l  +*a  +  l> 


(3-35) 


where  x„  is  the  "bridge"  point,  x,.,  the  end  point,  and  x,H ,  the 
second  point  at  the  beginning  of  the  record. 

DC  offset  and  rise  rate  errors  are  removed  next.  The  DC 
offset  is  a  function  of  the  amplitude  of  a  sampled  data  point 
and  it  is  independent  of  time  and  location.  The  amplitude  is 
scaled  by  a  constant  to  remove  this  error.  The  rise-rate 
error,  due  to  the  feedback  capacitor  charging,  is  time 
dependent.  It  is  zero  for  the  first  sample,  and  it  is  a 
linear  function  of  time  for  the  remaining  samples.  Once  the 
rise  slope  is  known,  the  error  can  be  subtracted  from  each 
sample  value  based  upon  the  length  of  time  it  occurs  relative 
to  the  first  data  point  in  the  record. 

Circuit  Layout 

The  output  from  the  charge  amplifier  has  level  gain  below 
1  Hz  and  has  a  calculated  3  dB  point  of  17  Hz.  This  feature 
implies  that  low  frequencies  will  pass  relatively  unaffected, 
but  that  high  frequencies  will  be  significantly  attenuated. 
The  film's  output  for  a  2  psi  load  will,  at  best,  be  0.3  V, 
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based  upon  a  110  pm  thick  film.  Thus,  the  charge  amplifier's 
output  signal  will  require  amplification  in  order  to  obtain 
adequate  resolution  with  the  A/D  board. 

An  instrumentation  amplifier  with  a  programmable  gain 
provides  accurate  gains  of  1,  10,  100,  and  1000,  and  it  can  be 
used  to  amplify  the  charge  amplifier’s  output.  An  Analog 
Device  AD524BD  instrumentation  amplifier  was  selected  to 
satisfy  this  requirement.  Figure  3-11  shows  the  arrangement 
of  the  charge  amplifier  with  the  instrumentation  amplifier. 

A  Qua  Tech  A/D  board  in  a  Zenith  Z-248  computer  was  used 
to  sample  and  record  the  data.  The  circuit  block  diagram  in 
Figure  3-12  shows  the  general  arrangement  used  for  testing  and 
deve 1 opment . 

PVDF  Transducer 

One  centimeter  square  pieces  of  film  were  cut  from  bulk 
pieces  of  PVDF  film  using  an  X-acto  knife.  The  bulk  film 
sheet  was  placed  on  a  large  piece  of  glass  to  minimize  damage 
to  the  electrode  metal  which  could  short  circuit  the  film  and 
make  it  unusable.  The  film  was  tested  for  electrical 
continuity  to  ensure  that  all  the  cuts  were  clean.  Next, 
shielded  two-wire  leads  were  attached  to  the  film  using 
h i gh 1 y-conduct i ve  silver  paint.  A  layer  of  adhesive  tape  was 
used  to  protect  the  film  from  corrosion  and  to  reinforce  the 
lead  attachments.  After  the  transducer  was  finished,  it  was 
bonded  to  a  cylinder  with  an  adhesive  spray  and  covered  with 
another  layer  of  adhesive  tape  to  prevent  debonding. 
Throughout  the  process,  the  film  was  not  touched,  but  handled 
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Figure  3-11.  Charge  Amplifier  with  Instrumentation  Amplifier 
Gain  Set  to  1000. 
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Test i ng . 


I 

I 
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with  tweezers  to  avoid  fingerprints  which  can  corrode  the 
electrode  metallization  (5:49). 

Test  and  Evaluation 

The  final  design  area  involved  the  development  of  a 
procedure  for  testing  the  circuit  and  film  transducer.  Each 
procedure  represented  a  small  step,  so  that  the  insight  gained 
would  contribute  to  the  next. 

Film  Characterization.  The  film  was  characterized  to 
determine  its  time  constant  and  sensitivity  S, .  That  is, 

1)  The  film  transducer  was  placed  between  two  pads  that 
matched  the  electrode  area  and  did  not  contact  the  lead 
attachment  area. 

2)  The  upper  pad  had  a  weight  placed  on  it  to  act  as  an 
inertial  mass  to  avoid  transients. 

3)  A  weight  of  known  mass  was  placed  on  the  inertial 
mass.  It  was  quickly  removed  by  hand  to  create  an 
output  proportional  to  the  mass.  As  the  mass  was 
removed,  the  output  voltage  would  rise,  and  the  film’s 
sensitivity  could  be  found  from  the  magnitude  of  the 
response  and  the  time  constant  from  the  charge  decay. 

4)  The  film's  output  was  fed  into  a  simple  voltage 
follower  shown  in  Figure  3-13.  The  same  operational 
amplifier  used  in  this  circuit  was  also  used  in  the 
charge  amplifier  network  so  that  the  input  impedances 
would  be  the  same  to  ensure  similar  response 
character i st i cs . 
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Figure  3-13.  Unity  Gain  Voyage  Follower 
(25:2-67)  . 


5)  The  film's  output  was  recorded  on  a  digitizing 
oscilloscope  (Tektronx  2430A)  set  in  its  s i ng 1 e- tr i gger 
mode.  Sample  responses  were  recorded  through  the  GP 1 B 
bus  for  subsequent  graphing  and  analysis. 

6)  This  procedure  was  repeated  27  times  with  three 
different  weights  (50,  100,  205  gm)  for  each  thickness 
of  film.  This  effort  facilitated  calculating  a 
statistical  average  for  each  weight  and  film 

comb i nat i on . 

Circuit  Character i zat ion .  There  are  four  required 
parts  to  characterize  the  charge  amplifier  circuit:  1)  rise 
rate  minimization;  2)  frequency  response;  3)  reset  pulse 
minimization;  and  4)  DC  offset  character i zat i on .  The  reset 
pulse  width  and  DC  offset  testing  are  i nterre 1 ated ,  but  they 
were  tested  individually. 
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Frequency  Response.  The  charge  amplifier's  frequency 
response  was  tested  with  the  reset  pulse  command  disabled  so 
that  complete  signals  could  be  measured.  That  is, 

1)  A  signal  generator  (Tektronix  FQ502)  was  connected  to 
the  input  of  the  charge  amplifier. 

2)  The  output  was  not  amplified  but  fed  directly  into  a 
digitizing  oscilloscope  (Tektronix  2430A). 

3)  A  sine  wave  of  known  frequency  and  amplitude  (measured 
on  the  oscilloscope)  was  connected  to  the  charge 
amplifier's  input,  and  the  circuit  output  was  measured. 

4)  The  frequency  was  swept  from  0.1  Hz  to  1  MHz.  The 
results  were  recorded  and  plotted  in  Figure  4-2  to 
produce  the  frequency  response  curve  for  a  resistive 
input . 

Rise  Rate  Character i zat i on .  The  circuit  rise  rate  defines 
the  linear  charge  rate  of  the  feedback  capacitor.  It  varies 
with  different  Cfk  and  Rfl  values.  The  rise  rate  creates  a 
linear  offset  error  which  starts  after  the  first  element  in 
the  data  record.  It  can  be  subtracted  during  the 

r econstruct i on  process  once  it  is  known.  That  is, 

1)  The  charge  amplifier  and  reset  inputs  were 
d i sconnected . 

2)  A  digitizing  oscilloscope  (Tektronix  2430A)  was 
connected  to  the  amplifier's  output. 

3)  The  voltage  supply  was  established  at  ±10V  DC. 

4)  The  feedback  capacitor  was  shorted  across  its 
terminals  to  discharge  it. 
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5)  The  output  was  recorded  on  the  oscilloscope  in  units 
of  volts  per  second. 

6)  The  feedback  resistor  and  capacitor  were  varied  to 
see  the  effect  on  rise  rate  by  repeating  steps  4  and 
5.  The  feedback  time  constant  (Rfk*Cfk)  target  was 
0.01  seconds. 

Pulse  Width  Minimization.  The  ideal  pulse  width  would  be 
very  short,  and  it  would  completely  discharge  the  feedback 
capacitor.  In  practice,  its  limit  is  determined  by  the  reset 
switch  resistance  and  feedback  capacitance.  As  the  reset  time 
is  reduced,  the  capacitor  does  not  fully  discharge,  and  a  DC 
offset  results.  If  the  reset  time  is  increased  beyond  5  time 
constants,  no  benefits  are  gained  since  the  capacitor  is  over 
99%  discharged.  That  is, 

1)  The  input  of  the  charge  amplifier  with  the  reset 
switch  was  connected  to  a  signal  generator  (Tektronix 
FG502),  and  a  pulse  generator  (Wave  Tek  278)  was 
connected  to  the  reset  circuit. 

2)  The  charge  amplifier's  output  was  connected  to  a 
digitizing  oscilloscope  (Tektronix  2430A). 

3)  The  pulse  amplitude  was  set  to  ±  10  V.  A  50  Hz,  ±  IV 
square  wave  was  connected  to  the  input  of  the 

amp  1 i f i er . 

4)  The  pulse  width  was  varied,  and  the  amplifier's 
output  was  recorded  on  the  oscilloscope. 

5)  The  output  was  triggered  on  the  reset  pulse,  and  the 
offset  was  measured  and  recorded. 
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DC  Offset  Character i zat ion .  The  DC  offset  values  results 


from  the  charge  decay  of  the  feedback  capacitor  and  charge 
injection  from  the  gate  to  the  sampling  mode  through  the  gate 
to  channel  capacitance  (14:100).  It  causes  a  DC  shift  in  the 
output  which  results  in  an  error  in  the  reconstructed  time 
history.  If  the  DC  shift  can  be  character  i  zed ,  it  can  be 
removed . 

This  test  used  the  same  procedure  as  the  pulse  width  test 
with  two  exceptions: 

1)  The  pulse  width  was  kept  constant  at  five  time 
constants  (40  ps);  and 

2)  The  square  wave  amplitude  was  varied  from  7  mv  peak- 
to-peak  to  1.75  V  peak-to-peak .  After  the  data  was 
collected,  a  statistical  average  of  the  maximum  and 
minimum  values  for  the  sampling  period  was  collected. 
The  values  were  plotted  to  obtain  the  DC  offset 
character izat ion . 

Secondary  cjain  Char acter  i  zat  i on  .  The  secondary  gain 
provides  the  amplification  needed  to  obtain  adequate 
resolution  on  the  A/D  board.  It  should  be  linear  at  all  gains 
below  32,000  Hz  to  ensure  accurate  results. 

1)  The  secondary  amplifier  (AD524BD)  input  was  zeroed  by 
grounding  the  inputs  and  setting  the  gain  to  1000. 

The  10  kohm  resistance  offset  nulling  potentiometer 
was  adjusted  so  that  the  output  was  zero  on  the 
osci 1 loscope . 
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2)  The  AD524BD's  output  was  zeroed  using  the  same 
procedure  described  in  step  1,  but  the  gain  was  set 
to  one,  and  the  output  potentiometer  was  adjusted  to 
zero  at  the  output. 

3)  A  signal  generator  (Tektronix  FG502)  was  connected 
to  the  inputs.  The  frequency  was  swept  from  1  to  1C 
MHz  with  a  sine  wave  of  constant  amplitude.  The  gain 
was  varied  from  1,  10,  100,  to  1000,  and  the 
frequency  sweep  was  repeated  each  time. 

4)  The  output  and  input  were  monitored  on  a  dual  trace 
oscilloscope.  The  data  were  recorded  and  reduced. 

Time  History  Reconstruction  Test 

In  order  for  the  piezoelectric  film  transducer  to  be 
useful,  an  accurate  time  history  had  to  be  r econstructed . 
This  test  would  determine  if  the  process  would  work.  That  is, 

1)  The  charge  amplifier  and  reset  circuit  was  configured 
with  a  500  Hz  reset  pulse  that  was  40  ps  wide. 

2)  The  input  was  connected  to  the  signal  generator 
(Tektronix  FG502)  producing  a  50  Hz  sine  wave  with  a 
known  amplitude. 

3)  The  output  of  the  signal  generator  and  charge 
amplifier  were  recorded  on  the  Qua  Tech  A/D  board 
at  a  rate  of  5000  samples  per  second. 

4)  The  data  were  post- processed  to  remove  the  DC  offset, 
rise  rate  error,  and  to  integrate  the  samples.  The 
results  were  plotted  with  the  original  input,  the 
"chopped"  output,  and  the  reconstructed  output. 
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System  Testing 


This  test  would  determine  if  a  piezoelectric  film  pressure 
transducer  is  a  viable  option  for  measuring  pressure 
frequencies  ranging  from  DC  to  2,500  Hz.  The  film  response 
was  compared  to  the  output  of  a  calibrated  dynamic  pressure 
transducer  (Endevco  8510-B).  Each  film  thickness  was  tested 
as  a  separate  run  due  to  the  2  channel  limit  of  the  A/D  board. 
That  is, 

1)  The  three  piezoelectric  traducers  were  mounted  along 
the  side  of  a  dynamic  pressure  gage  as  shown  in 
Figure  3-14.  It  was  assumed  that  the  lateral 
distance  from  the  "reference"  gage  would  not  create 

a  significant  pressure  difference  in  a  large  volume 
flow  field. 

2)  The  output  of  the  reference  gage  was  connected  to  its 
signal  conditioner. 

3)  The  charge  amplifier  was  driven  with  a  500  Hz  reset 
rate  with  a  100  ps  pulse  width. 

4)  The  outputs  of  the  transducers  were  zeroed,  and  a 
nominal  condition  was  recorded  on  the  A/D  board  at 
5000  samples  a  second.  This  would  facilitate  an 
additional  offset  to  be  removed  from  the  test  data. 

5)  The  cylinder  with  the  gages  was  placed  in  a  flow 
field  and  the  flow  around  the  cylinder  was  measured 
with  a  pitot  tube.  This  measurement  established  the 
relative  magnitude  of  the  flow  field.  Figure  3-15 
shows  the  cylinder  and  flow  field  arrangement. 
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Figure  3-15.  Cylinder  with  Pressure  Transducers  in  the  Flow 


Field. 


I 

I 
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6)  A  test  run  was  initiated  by  blocking  the  flow  with 

a  solid  aperture.  The  data  acquisition  process  was 
initiated  and  the  plate  was  removed.  This  condition 
created  a  "no-flow"  and  a  "ramp"  pressure 
distribution  on  the  transducers.  The  "no-flow" 
condition  is  important  since  a  nominal  condition  must 
be  known  in  order  to  have  an  accurate  starting  point 
for  the  time  history  r econstruct i on  process. 

7)  The  data  from  the  "reference"  gage  was  scaled  with 
its  calibration  data.  The  piezoelectric  gage  data 
was  reconstructed  with  the  same  process  used  in  the 
reconstruct  ion  test. 

8)  The  results  for  each  gage  were  put  on  the  same  plot 
and  analyzed. 

Summary 

This  chapter  presented  additional  material  required  to 
design  and  understand  the  charge  amplifier  reset  circuit. 
Design  calculations  were  presented  to  determine  an  initial 
point  for  the  performance  testing  and  development.  Finally, 
the  test  procedure  was  presented  to  verify  and  improve  the 
transducer's  design. 
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IV.  Experimental  Data  and  Analysis 


I ntroduct i on 

The  presentation  of  this  section  parallels  the  order  of 
the  tests.  The  results,  conclusions,  and  additional  insights 
gained  from  testing  are  discussed. 

Film  Character i zat ion 

The  results  from  the  sensitivity  tests  are  presented  in 
Table  4-1.  The  PVDF  film's  response  in  all  cases  was  within 
7%  of  the  theoretical  values.  The  relatively  large  standard 
deviations  can  be  attributed  to  the  experimental  test 
variations.  If  the  weight  was  not  removed  at  the  same  speed 
or  perfectly  vertical ,  a  si ight  deviation  from  the  ideal  value 
was  noted.  The  standard  deviation  was  found  from  the 
following  relationship: 


Each  measurement  was  repeated  27  times  to  establisn  a  simple 
average . 

The  measured  charge  decay  rate,  or  film  time  constant  was 
significantly  different  compared  to  the  calculated  result. 
The  theoretical  •  alue  was  on  the  order  of  30  minutes;  whereas, 
the  actual  values  were  on  the  order  of  1.5  seconds.  Also,  the 
longest  decay  should  have  been  associated  with  the  28  pm  thick 
film,  but  the  observed  decay  was  longest  for  the  110  pm  thick 
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Table  4-1.  PVDF  Film  Sensitivity. 


F i lm-Property 

50  gm 

100  gm 

206  gm 

28  pm-S,  Theory 

28.5  mv 

57.0  mv 

117.0  mv 

28  pm-S,  Measured 

30.4  mv 

59.3  mv 

122.3  mv 

28  pm-Variance 

5.9  mv 

8.5  mv 

17.3  mv 

52  pm-S,  Theory 

52.9  mv 

105.6  mv 

217.0  mv 

52  pm-S,  Measured 

51.3  mv 

107.0  mv 

212.5  mv 

52  pm-Variance 

6.9  mv 

10.9  mv 

17.6  mv 

110  pm-S,  Theory 

116.6  mv 

223.4  mv 

458.4  mv 

110  pm-S,  Measured 

119.4  mv 

220.0  mv 

463.0  mv 

110  pm-Variance 

12.8  mv 

18.5  mv 

27.5  mv 

film.  Figure  4-1  shows  a  sample  charge  decay  plot  for  the 
110  pm  film. 

The  difference  between  the  observed  and  theoretical  charge 
decay  times  can  be  attributed  to  non- ideal  components  and  the 
small  amount  of  charge  generated.  Anytime  an  electrical 
interface  with  a  finite  impedance  is  connected  to  a  capacitor, 
it  will  sink  the  storage  charge,  thus  increasing  the  charge 
decay  rate  associated  with  of  the  capacitor.  One  method  of 
more  accurately  measuring  the  charge  decay  rate  would  be  to 
use  a  higher  input  impedance  instrument  such  as  an 
electrometer  (which  was  not  available). 

The  PVDF  film  characterization  also  revealed  the  polarity 
of  the  film.  This  would  allow  the  transducers  to  be 
constructed  so  that  a  compressive  force  creates  a  positive 
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voltage  deflection.  This  is  the  convention  commonly  used  by 
the  "reference"  pressure  transducer. 

Circuit  Character i zat ion 

This  set  of  tests  included  the  evaluation  of  the  charge 
amplifier  frequency  response,  rise  rate  character i st i c s , 
pulse-width  minimization,  and  DC-offset  character i st i cs . 

Frequency  Response.  The  gain  of  the  charge  amplifier  was 
plotted  in  dB  [  20*LOG,,(  Vu,/Vj, )  ]  as  a  function  of  frequency  in 
Figure  4-2.  The  calculated  3  dB  point  was  17  Hz,  and  the 
measured  result  occurred  at  32.5  Hz.  The  gain  decreased 
monoton i ca 1 1 y  as  the  frequency  increased.  The  linear  least- 
squares  line  reveals  that  the  relationship  between  gain  and 
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Figure  4-2.  Frequency  Response  of  the  Charge  Amplifier 
Circuit. 


frequency  is: 


Gain  =  21.77  -  7.63  *  log  (Frequency) 


(4-1) 


The  behavior  of  the  linear  response  establishes  excellent 
agreement  between  the  theory  and  the  experimental  results. 
The  rapid  fall-off  feature  at  low  frequencies  affords  the 
circuit  its  integration  properties  for  resistive  inputs. 
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Rise  Rate  Character i st i cs .  The  rise  rate  measurements 
revealed  that,  as  the  size  of  Cfk  was  decreased,  the  charging 
rate  increased.  These  results  are  summarized  in  Table  4-2. 
Using  a  larger  Rfk  increased  the  rise  rate  slightly.  Different 
combinations  of  R  and  C  with  the  same  RC  product  resulted  in 
different  rise  rates.  In  each  case,  Cfk  was  the  dominant 
parameter.  To  keep  the  RC  product  at  approximately  0.01  (1/s) 
and  to  minimize  rise  rate  error,  C(k  was  set  at  0.02  pF . 
Consequently,  Rfl  was  established  at  470,000  ohms  to  improve 
gain,  since  gain  is  proportional  to  1/RC. 


Table  4-2.  Rise  Rates  for  the  Circuit's  Various  RC 

Combinations . 


cfk  (nF) 

Rri  (kQ) 

Rise  Rate  (V/S) 

RC  (1/s) 

0.002 

10 

0.46 

2E-5 

0.002 

100 

! 

0.53 

2E-4 

0.002 

1000 

0.55 

2E-3 

0.01 

10 

0.112 

IE-4 

0.01 

1000 

0.120 

IE-2 

0.02 

10 

0.052 

2E-4 

0.02 

470 

0.057 

9.4E-3 

0.02 

1000 

0.062 

2E-2 

Pulse  Width  Minimization.  The  pulse-width  variations 
investigated  were  found  to  increase  the  DC  offset  when  it  was 
set  to  less  than  5  time  constants.  In  fact,  if  the  pulse  time 
constant  was  measured,  it  was  found  to  be  within  1  ps  of  the 
calculated  value.  Figure  4-3  shows  a  reset  pulse  with  the 
associated  switching  noise  and  feedback  capacitor  charge 
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decay.  As  the  pulse-width  was  increased  to  be  larger  than  5 
time  constants,  there  was  no  additional  decrease  in  the  DC 
offset.  As  a  result,  the  reset  pulse-width  was  selected  to  be 
equal  to  5  time  constants. 

DC  Offset  Character i st i cs .  The  magnitude  of  the  last 
reset  period  sample  relative  to  the  first  reset  period  sample 
(DC  offset)  was  found  to  vary  linearly  with  magnitude.  This 
means  that  the  DC  offset  could  easily  be  removed  during  the 
time  history  data  reconstruct  ion  process.  Figure  4-4  shows  a 
data  sample  used  to  determine  a  DC  offset  point.  The  complete 
results  are  plotted  in  Figure  4-5.  The  linear  least-squares 
line  has  the  equation: 
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<r 


Y-  0 . 04904X  +  6 . 58i?-4 


(4-2) 


Secondary  Gain.  The  AD524BD  instrumentation  amplifier  is 
designed  to  give  precise  gains  of  1,  10,  100,  and  1000.  The 
test  results  agreed  precisely  with  the  specification  values 
throughout  the  entire  operational  frequency  range;  Table  4-3 
shows  the  cutoff  frequency  for  each  level  of  gain. 

Table  4-3.  Measured  and  Specified  Bandwidth  for  AD524BD 

Instrumentation  Amplifier. 


Qain 

Measured  -3  dB  Point 

Specified  -3  dB  Point  j 

1 

376  kHz 

1  MHz 

10 

430  kHz 

400  kHz 

100 

309  kHz 

150  kHz 

1000 

43  kHz 

25  kHz 

Figure  4-4.  Sample  DC  Offset  Data. 
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Figure  4-5.  DC  Offset  as  a  Function  of  Input  Voltage 


Time  History  Reconstruction  Process 

With  the  DC  offset  characterized,  the  "chopped"  time  history 
could  be  reconstructed.  Figure  4-6  shows  a  plot  of  a 
"chopped"  sine  wave,  the  reconstructed  sine  wave,  and  the 
input  sine  wave.  The  original  chopped  sine  wave  had  a  nonzero 
mean.  To  facilitate  the  r econstruct ion  process,  the  mean  was 
adjusted  to  zero  by  subtracting  the  operational  amplifier's 
induced  offset  from  all  data  points  (not  the  DC  offset  value). 
The  nonzero  mean  was  caused  by  an  imprecise  trim  adjustment  of 
the  charge  amplifier. 

In  this  simple  case,  the  reset  points  were  located  using 
a  saturated  condition  metric;  that  is,  by  counting  10  points 
to  the  next  reset  period,  or  by  an  obvious  mismatch  in 
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Reconstructed  Sine  Wave 
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Figure  4-6.  Reconstructed  Sine  Wave  Showing  the  Sampled  Data 
and  Input  Sine  Wave. 


amplitude.  When  the  samples  were  integrated,  the  "bridge” 
points  showed  the  s 1 i ght  mi smatch .  Various  other  methods  were 
used  to  try  to  make  the  transition  smoother.  In  some  cases 
the  "bridge"  was  smooth,  but  at  other  points,  the  "bridge"  was 
worse.  Therefore,  the  point  averaging  method  appeared  to  be 
the  most  effective  procedure. 

The  DC  offset  reductions  significantly  decreased  a 
steadily  increasing  amplitude.  However,  the  amplitude  still 
was  not  level.  It  still  manifested  an  amplitude  rise  due  to 
the  circuit's  rise  rate.  The  DC  offset  adjusting  factor  was 
increased  and  decreased  to  see  if  it  would  correct  the 
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amplitude  rise.  This  consideration  contributed  no 
improvements,  so  a  linear  least-squares  line  was  fitted 
through  the  peak  amplitude  points,  and  the  circuit's  rise  rate 
slope  was  determined,  The  slope  (-1.172673  volts  per  second) 
was  subtracted  from  each  point  as  a  function  of  its  time  shift 
relative  to  the  first  recorded  point.  This  data  processing 
procedure  made  the  amplitude  consistent. 

The  reconstructed  sine  wave  was  phase  shifted  by  nearly  90 
degrees  relative  to  the  input  signal.  Since  the  signal  was  50 
Hz,  it  was  attenuated  by  10  dB,  as  predicted  by  Figure  4-2. 
This  demonstration  shows  that  the  reconstruction  technique  was 
accurate  for  the  test  case  with  a  known  input. 

System  Testing 

The  laboratory  instrumentation  sampling  rate  of  5000  Hz 
limits  the  frequency  content  analysis  to  2500  Hz  due  to  the 
Nyquist  sampling  theorem.  With  this  limitation,  the  data  were 
digitally  filtered  with  an  8-pole,  0.01  dB  ripple,  low-pass 
Chebyshev  f i Her  to  prevent  aliasing  in  the  frequency  domain. 
The  filter  not  only  reduce  aliasing,  but  also  alternates  the 
20  kHz  signal  that  was  present  in  all  the  data  when  viewed  on 
the  oscilloscope.  The  source  of  the  20  kHz  signal  source 
could  not  be  found,  so  it  was  filtered. 

There  was  also  a  60  Hz  noise  signal  in  all  of  the  measured 
data  due  to  a  ground  loop.  The  ground  loop  occurred  between 
the  pulse  generator  and  the  circuit  ground.  Whenever  the 
pulse  generator  was  completely  disconnected,  the  60  Hz  signal 
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disappeared.  In  some  cases,  the  60  Hz  noise  was  as  large  as 
40  mv  RMS.  Various  grounding  techniques  were  implemented,  but 
none  contributed  an  improvement.  Consequently,  the  dynamic 
pressure  transducer  data  was  filtared  with  a  60  Hz  4-pole, 
-40  dB  stopband  filter  before  any  analysis  or  data  processing 
was  accomplished.  The  sampled  data  was  filtered  a'ter 
r econstruct i on ,  but  before  additional  processing  was 

imp  1 emented . 

Frequency  Content.  Power  spectral  dersity  plots  were 
calculated  for  both  the  PVFD  film  data  and  the  pressure  gage 
data  to  determine  their  frequency  content.  This  behavior 
indicates  that  the  PVDF  film  was  responding  dynamically  and 
similar  to  the  commercial  pressure  gage.  The  frequency 
resolution  was  on  the  order  of  20  Hz,  so  accurate  low 
frequency  analysis  was  not  possible  from  the  PSDs. 
Consequently,  the  low  frequency  magnitude  response  was 
determined  from  the  time  history. 

Comparing  the  two  PSDs  shows  that  the  film's  frequency 
response  matched  the  gage  response  above  200  Hz.  Below  200 
Hz,  the  film  PSD  shows  that  the  response  is  40  dB  greater  than 
the  relative  magnitude  of  the  gage  PSD.  Figures  4-7  and  4-8 
illustrate  the  two  representative  PSDs  fo>  the  110  pm  thick 
film.  Between  200  and  1000  Hz,  the  peaks  of  the  PVDF  film’s 
PSD  occur  close  to  the  same  frequencies  as  those  of  the  gage. 
At  700  Hz,  the 

film  response  is  smaller  relative  to  the  gage,  and  it  is  20 
dB  less  compared  to  the  other  peaks  on  the  film’s  PSD. 
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Power  (votts’Vota/Hz) 


Figure  4-7.  PSD  of  Dynamic  Pressure  Gage  Time  History  Data. 
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Between  1000  and  2500  Hz,  the  responses  are  not  nearly  as 
similar.  The  gage  response  decreases  by  more  than  20  dB  after 
1000  Hz,  and  then  it  increases  40  dB  and  peaks  at  2060  Hz,  and 
then  decreases  after  this.  The  film  shows  a  strong  response 
up  to  1160  Hz,  and  then  it  has  a  20  dB  decreases  out  to  1200 
Hz.  Above  1200  Hz,  the  PVDF  film's  response  increases  by  20 
dB  and  peaks  at  2200  Hz.  Both  responses  decrease  by  more  than 
220  dB  after  2500  Hz  due  to  filtering. 

Below  200  Hz  the  two  responses  differ  considerably.  The 
PVDF  film's  response  shows  a  dominant  low  frequency  component 
at  2  Hz  that  diminishes  as  the  frequency  increases  to  180  Hz. 
The  gage  frequency  response  decreases  steadily  from  2  to  60 
Hz,  the  stopband’s  center  frequency,  and  it  increases  by 
almost  60  dB  and  peaks  at  160  Hz.  The  magnitude  of  the  peak 
on  the  gage  response  is  the  significant  difference  between  the 
film  and  gage.  The  PVDF ' s  film  response  manifests  a  stronger 
low  frequency  response,  and  the  gage  indicates  a  large 
response  centered  at  160  Hz. 

Time  History.  The  time  histories  of  the  film  and  gage 
show  similar  trends.  These  are  shown  in  Figures  4-9  and  4-10. 
Overall,  the  gage  shows  a  much  greater  dynamic  response,  but 
the  low  frequency  and  DC  trends  of  both  gages  are  very 
similar.  On  a  point-by-point  basis,  the  film's  response 
appears  to  be  time  shifted  by  approximately  5  ms.  This  time 
shift  is  attributed  to  be  the  cause  of  the  phase  shift 
observed  in  the  charge  amplifier's  output. 
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Figure  4-9.  Filtered  Data  for  the  Dynamic  Pressure  Gage 


Figure  4-10.  Filtered  and  Corrected  PVDF  Film  Gage  Data. 
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Figures  4-11  and  4-12  show  the  time  histories  for  52 
thick  film.  The  PVDF  film’s  response  clearly  shows  the 
pressure  rise  due  to  the  unblocking  of  the  air  flow.  This 
indicates  that  the  PVDF  film  transducer  can  detect  low 
frequency  pressure  fluxuations.  The  dynamic  pressure  gage 
does  not  show  a  similar  step  response  because  it  does  not  give 
an  integrated  output. 


Filtered  Dynamic  Pressure  Gage  Data 
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Figure  4-11.  Filtered  Dynamic  Pressure  Gage  Data. 
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Figure  4-12.  Filtered  and  Corrected  52  Micron  Thick  PVDF  Film 
Gage  Data. 

Observations 

Several  areas  were  not  evaluated  such  as  pyroelectric 
effects,  electromagnetic  inter f erence ,  and  external 
vibrational  effects.  In  the  PVDF  film’s  response  tests,  if  a 
warm  hand  was  brought  near  any  film  sample,  it  would  produce 
several  hundred  millivolts,  which  is  as  large  as  the  response 
to  a  pressure.  Thus,  the  pyroelectric  effect  requires 
additional  investigation.  The  effect  of  electromagnetic 
radiation  was  not  apparent  in  any  testing.  However,  external 
vibration  sources  can  be  a  significant  problem  with  the 
piezoelectric  film.  If  something  was  dropped  on  the  test 
bench,  the  PVDF  film,  but  not  the  commercial  gage,  would  show 
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a  response.  It  would  also  respond  to  any  stresses  or 
vibrations  in  the  model  or  surface  it  was  mounted  on.  The 
next  section  will  address  some  possible  improvements  for  these 
prob 1 ems . 

Cone  1  us i on 

All  of  the  performance  tests  proved  to  be  beneficial  and 
facilitated  the  design  process  and  the  choice  of  components. 
The  time  history  r econstruct i on  test  proved  that  an  accurate 
response  picture  can  be  obtained  from  the  sampled  data.  This 
process  facilitated  collecting  a  pressure  history  of  the  film 
and  its  reconstruct  ion  to  determine  if  the  PVDF  pressure 
transducer  was  a  viable  concept.  The  initial  data  analysis 
shows  agreement  between  the  film  and  gage  responses. 
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V. 


Conclusions  and  Recommendations 


Cone  1  us i ons 

The  purpose  of  this  investigation  was  to  design  and 
evaluate  the  performance  of  a  piezoelectric  film  pressure 
transducer  to  measure  frequencies  from  DC  to  over  32  kHz.  To 
do  this,  a  reset  and  sampling  technique  was  implemented  with 
an  analog  circuit  design.  A  process  had  to  be  developed  and 
verified  to  reconstruct  the  sampled  segments  of  data  into  a 
complete  record.  Next,  the  circuit  and  reconstruct i on  process 
was  used  to  determine  if  piezoelectric  film  could  be  applied 
to  a  body  to  function  as  a  pressure  transducer. 

This  investigation  was  successful  concerning  the  design 
and  performance  evaluation  of  the  analog  reset  circuit  and 
sampled  data  reconstruction  process  for  a  known  input.  The 
piezoelectric  film  gage  appears  to  have  potential  as  a 
pressure  gage,  but  additional  development  is  required  before 
it  could  be  used  in  a  practical  application.  The  film  agreed 
with  the  dynamic  pressure  gage  at  frequencies  below  1000  Hz. 
The  dynamic  response  of  the  piezoelectric  film  was  limited 
compared  to  the  dynamic  pressure  gage  at  frequencies  above 
1000  Hz. 

Recommendations 

The  following  recommendations  are  presented  to  improve  the 
design  and  development  of  the  piezoelectric  film  pressure 
transducer . 
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1)  The  charge  amplifier  circuit  needs  to  be  better  isolated 
from  the  60  Hz  grounding  loop.  By  making  a  differential 
amplifier  with  several  operational  amplifiers,  the  film 

and  amplifier  would  not  be  connected  to  the  same  ground  as 

the  operating  instrumentation.  This  feature  would  facilitate 

avoiding  ground  loops  and  minimize  the  effect  of  60  Hz  noise. 

2)  With  a  more  sophisticated  A/D  device,  the  reset  command 
pulse  and  data  sampling  process  could  be  synchron i zed . 
This  would  eliminate  ambiguity  in  identifying  the  start  of 
each  sampling  period.  As  a  result,  the  first  data  point 
of  each  sample  could  be  used  and  this  feature  would 
improve  accuracy.  In  addition,  the  data  reduction  process 
would  be  simpler  and  faster. 

3)  The  effects  of  vibration  and  heat  on  the  piezoelectric 
film  need  to  be  investigated.  In  a  more  advanced  test 
configuration,  a  "reference"  gage  that  was  isolated  from 
the  pressure  field,  could  be  used.  The  signal  of  the 
reference  gage  could  be  subtracted  from  the  pressure  gage, 
thus  nulling  any  p yro,  vibration,  or  strain  effects. 

4)  A  method  must  be  developed  to  identify  the  sensor's 
initial  condition  on  a  long  data  run.  Possibly,  a  slow 
sampling  and  reset  rate  could  be  used  until  the  system 
stabilized,  then  the  rate  increased  to  get  the  desired 
data . 

5)  All  of  the  components  used  in  this  investigation  were  the 
best  available.  One  of  the  original  intents  was  to  use 
inexpensive  commonly  available  components.  Repeating  a 
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standard  test  with  inexpensive  components  would  determine 
the  effect  these  parts  have  concerning  charge  leakage  and 
accuracy.  As  a  result,  if  there  was  no  degradation  in 
performance,  these  parts  could  be  used. 

6)  Since  the  film  is  flexible  and  compliant,  it  can  be 
attached  to  almost  any  surface  with  an  adhesive,  for 
example,  an  area  on  a  aerodynamic  body  covered  with  a 
single  piece  of  film  with  an  etched  array  of  piezoel ectr  i  c 
gages  on  it.  This  arrangement  would  facilitate  the 
creation  of  a  pressure  map.  This  particular  area  needs  to 
be  fully  exploited  and  developed  to  capture  the  potential 
of  p i ezoe 1 ectr i c  film  as  a  pressure  transducer. 
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Table  A-1.  Materials 


Kynar  PVDF  film  (28  pin,  52  pm,  and  110  pm), 
Pennwalt  Corporation, 

Va 1  ley  Forge  PA 
AD549JH  operational  Amplifier, 

Analog  Devices, 

Norwood  MA 

AD524BD  Instrumentation  Amplifier, 

Analog  Devices, 

Norwood  MA 

0.02  pF  Teflon  Capacitor  ( J 1 1 B203KSW) , 
Component  Research  Company, 

Santa  Monica  CA 
3N170  Transistors 
1  MQ  Resistor 
470  k$2  Resistor 
10k  Q  Resistor 
0.1  pF  Ceramic  Capacitors 
10  kQ  Potentiometer 
1N914  Diode 

Conductive  Silver  Paint 
Spray  Adhesive 
2"  Polyester  Tape 
Twin  Lead  Shielded  Cable 
18  Gauge  Solid  Copper  Wire 
BNC  Test  Leads 


Table  A-2.  Equipment 


Tektronix  2430A  Digital  Oscilloscope, 

Tektronix,  Beaverton  OR 
Tektronix  FG502  Function  Generator, 

Tektronix,  Beaverton  OR 
Tektronix  DM501A  Multimeter, 

Tektronix,  Beaverton  OR 
Power  Ten  3103-A-2000  DC  Power  Supply, 

Power  Ten,  Sunnyvale  CA 
Wavetek  278  Function  Generator, 

Wavetek  San  Diego, 

San  Diego  CA 

Zenith  248  Data  Collection  Computer  (With  IEEE-488  Interface), 
Zenith  Data  Systems, 

Glenview  IL 

Qua  Tech  12  Bit  A/D  converter  System  (In  Z-248), 

Qua  Tech,  Akron  OH 
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Endevco  8510-B  Pressure  Transducer, 

Endevco , 

San  Juan  Capistrano  CA 
Endevco  4225  and  4423  Signal  Conditioners, 
Endevco , 

San  Juan  Capistrano  CA 
Weights  ( 3 ) 

Protoboar d 
X-acto  knife 
Tweezers 

Glass  Cutting  Surface 
2"  Acrylic  Tubing  (cylinder) 

100  psi  Regulated  Air  Supply 
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Appendix  B:  Recommended  Materials 


Table  B-1.  Recommended  Operational  Amplifiers. 

Analog  Devices  AD549 
Burr  -  Brown  OPA128 
PM  I  OP-80 


Table  B-2.  Recommended  Transistors. 


3N169 
3N170 
3N17  1 
Motoro 1  a 
Motoro  1  a 


2N7000 

VN2222LL 
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Appendix  C:  Data  Collection  and  Analysis  Programs 


Oscilloscope  with  IEEE-488  Interface  Data  Collection  Program 


Declarations  for  Quick  BASIC  4. 0/4. 5  applications. 

'  NOTE  :  include  this  file  only  if  you  are  usir.*  Quick  BASIC 
Version  4.0  or  higher. 

'  Common  GP I B  status  variables 

COMMON  SHARED  /NIST ATBLK/  IBSTA*,  IBERR*,  I BCNT* 

GP  I B  Subroutine  Declarations 


DECLARE 

SUB 

1  BBNA 

(BD*, 

BDNAME$ ) 

DECLARE 

SUB 

1  BCAC 

(BD*, 

V*) 

DECLARE 

SUB 

1  BCLR 

(BD*) 

DECLARE 

SUB 

1  BCMD 

(BD*, 

CMD$  ) 

DECLARE 

SUB 

1 BCMDA 

(BD* 

,  CMD$ ) 

DECLARE 

SUB 

1  BDMA 

(BD*, 

V*) 

DECLARE 

SUB 

1  BEOS 

(BD*, 

V*) 

DECLARE 

SUB 

1  BEOT 

(BD*, 

V*) 

DECLARE 

SUB 

IBFIND 

( 8DNAME$ ,  BD*) 

DECLARE 

SUB 

1  BGTS 

(BD*, 

V*) 

DECLARE 

SUB 

IBIST 

(BD*, 

V*) 

DECLARE 

SUB 

1  BLOC 

(BD*) 

DECLARE 

SUB 

1  BONL 

(BD*, 

V*) 

DECLARE 

SUB 

IBPAD 

(BD*, 

V*) 

DECLARE 

SUB 

1  BPCT 

(BD*) 

DECLARE 

SUB 

IBPPC 

(BD*, 

V*) 

DECLARE 

SUB 

IBRD  (BD*,  1 

RD$) 

DECLARE 

SUB 

1  BRDA 

(BD*, 

RD$  ) 

DECLARE 

SUB 

1  BRDF 

(BD*, 

FLNAME$ ) 

DECLARE 

SUB 

1  BRDI 

(BOSS, 

1 ARR* ( ) ,  CNT*) 

DECLARE 

SUB 

IBRDI A 

(BD* 

,  1 ARR*( ) ,  CNT*) 

DECLARE 

SUB 

1  BRPP 

(BD*, 

PPR*  ) 

DECLARE 

SUB 

IBRSC 

(BD*, 

V*) 

DECLARE 

SUB 

1  BRSP 

(BD%, 

SPR*) 

DECLARE 

SUB 

1  BRSV 

(BD*, 

V*) 

DECLARE 

SUB 

1  BSAD 

(BD*, 

V*) 

DECLARE 

SUB 

1  BS  1C 

(BD*) 

DECLARE 

SUB 

1  BSRE 

(BD%, 

V*) 

DECLARE 

SUB 

1 BSTOP 

(BD*) 

DECLARE 

SUB 

1  BTMO 

(BD*, 

V*) 

DECLARE 

SUB 

1  BTRAP 

(MASK*,  MODE*) 

DECLARE 

SUB 

1  BTRG 

(BD*) 

DECLARE 

SUB 

1 BWA 1 T 

(BD* 

,  MASK*) 

DECLARE 

SUB 

1  BWRT 

(BD*, 

WRT$) 

DECLARE 

SUB 

1  BWRTA 

(BD* 

,  WRT$) 
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DECLARE  SUB  I BWRTF  (BDSS,  FLNAME$) 

DECLARE  SUB  IBWRTI  ( BDSS ,  IARRSS(),  CNTSS) 
DECLARE  SUB  IBWRTIA  (BDSS,  IARRS5(),  CNTSS) 

'  GP I B  Function  Declarations 


DECLARE 

FUNCT 

ION 

1 LBNASS 

(BDSS, 

BDNAME$) 

DECLARE 

FUNCT 1  ON 

1 LCACSS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

1 LCLRSS 

(BDSS) 

DECLARE 

FUNCT 

ION 

1  LCMDSS 

(BDSS, 

CMD$,  CNTSS ) 

DECLARL 

FUNCT 

ION 

ILCMDASS 

(BDSS 

,  CMD$ ,  CNTSS) 

DECLARE 

FUNCT 

ION 

1 LDMASS 

(  BDSS , 

VSS) 

DECLARE 

FUNCT 

ION 

1  LEOSSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

ION 

1 LEOTSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

1  ON 

1  LF  1  NDSS 

(BDNAME$) 

DECLARE 

FUNCT 

ION 

1  LGTSSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

ION 

1 L 1 STSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

ION 

ILLOCSS 

(  BDSS) 

DECLARE 

FUNCT 

ION 

1  LONLSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

ION 

1  LPADSS 

C  BDSfe , 

VSS) 

DECLARE 

FUNCT 

ION 

1  LPCTSS 

(BD*) 

DECLARE 

FUNCT 

ION 

1  LPPCSS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

ILRD*  ( 

BDSS,  1 

RD$ ,  CNTSS) 

DECLARE 

FUNCT 

ION 

1  LRDASS 

(BDSS, 

RD$,  CNTSS) 

DECLARE 

FUNCT 

ION 

ILRDFSS 

(BDSS, 

FLNAME$) 

DECLARE 

FUNCT 

ION 

ILRDlSS 

(BDSS, 

1  ARRSS (  )  ,  CNTSS ) 

DECLARE 

FUNCT 

ION 

1  LRD  1  ASS 

(BOSS 

,  IARRSS(),  CNTSS) 

DECLARE 

FUNCT 

ION 

1  LRPPSS 

(BDSS, 

PPRSS) 

DECLARE 

FUNCT 

ION 

1 LRSCSS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

1  LRSPSS 

(BDSS, 

SPRSS ) 

DECLARE 

FUNCT 

ION 

1  LRSVSS 

(BDSS, 

VSS) 

DECLARE 

FUNCT 

ION 

1 LSADSS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

1 LS ICSS 

(BDSS) 

DECLARE 

FUNCT 

ION 

1  LSRESS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

1  LSTOPSS 

(BCS5) 

DECLARE 

FUNCT 

ION 

1 LTMOSS 

(BDSS, 

vss) 

DECLARE 

FUNCT 

ION 

ILTRAPSS 

(MASKSS,  MODESS) 

DECLARE 

FUNCT 

ION 

1  LTRGSS 

(BDSS) 

DECLARE 

FUNCT 

ION 

ILWAITSS 

(BDSS 

,  MASKSS) 

DECLARE 

FUNCT 

ION 

1 LWRTSS 

(BDSS, 

WRT$,  CNTSS) 

DECLARE 

FUNCT 

ION 

ILWRTASS 

(BDSS 

,  #RT$,  CNTSS) 

DECLARE 

FUNCT 

ION 

1 LWRTFSS 

(BDSS 

,  FLNAME$ ) 

DECLARE 

FUNCT 

ION 

1  LWRT  1  SS 

(BDSS 

,  IARRSS(),  CNTSS) 

DECLARE 

FUNCT 

ION 

1  LWRT  1  ASS  (BDSS,  IARRSS(),  CNTSS) 

INPUT  "NAME  DATA  FILE  FOR  DATA  STORAGE : XXXXX  NO  EXTENSION" 
NA$ 

NA$  =  "D:"  +  NA$  +  ".DAT" 

OPEN  NA$  FOR  OUTPUT  AS  #1 

UDNAME$  =  "DEVI" 

CALL  I  BF  I  ND('JDNAME$,  DEV  1  % ) 

WRT$  =  "CLE" 
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JU$  =  SPACE$( 50 ) 

CALL  I BWRT ( DEV  1  % ,  WRT$) 

CALL  I BRD ( DEV  1 % ,  JU$) 

WRT$  =  "CH 1 ?VOL " 

CALL  I BWRT ( DEV  1 % ,  WRT$) 

VO$  =  SP ACE$ ( 100) 

CALL  I BRD ( DEV  1  % ,  VO$) 

VO!  =  VAL ( VO$ ) 

VO!  =  VO!  *  .04 
PRINT  VO! 

WRT$  =  "HOR?BSE" 

CALL  I BWRT (DEV  1 % ,  WRT$) 
ho$  =  SPACES ( 100) 

CALL  I BRD ( DEV  1  % ,  ho$) 
ho  =  VAL(ho$) 
ho  =  (ho  *  20)  /  1024 
PRINT  ho 

WRT$  =  "PAT  OFF" 

CALL  I BWRT ( DEV  1% ,  WRT$) 

INPUT;  "PRESS  ENTER  TO  COLLECT  WAVEFORM  "  ;  XXX$  'WRT$="  WA?" 
WRT$  =  "CURV? " 

CALL  I  BWRT ( DEV  1$ ,  WRT$) 

RD$  =  SPACES ( 5100) 


CALL  I BRD ( DEV  1% ,  RD$ ) 

HOI  =  0 

WRT$  =  "LOC  OFF" 

CALL  I  BWRT  (DEV  196 ,  WRT$) 

Z  =  0 
K  =  0 

txtlength  =  LEN ( RD$ ) 

p  =  0 

FOR  I  =  1  TO  txtlength 

C$  =  MID$(RD$,  I  ,  1 ) 

IF  C$  =  THEN 

Z  =  VAL(D$):  D$  =  C$  =  PRINT  #1,  USING 

1"#.#### . ;  HOI; 

PRINT  #1,  Z  *  VO!:  HOI  =  HOI  +  ho 

p  =  p  +  1 

END  IF 

D$  =  D$  +  C$ 

K  =  K  +  1 

NEXT  ' 

PRINT  p,  K 
CLOSE 

This  program  was  written  in  Microsoft  Quick  Basic  version  4.5. 
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Qua  Tech  A/D  Board  Data  Collection  Program 


10  CL  S 
20  COLOR 
30  PRINT 
40  PRINT 
50  PRINT 
60  PRINT 
70  PRINT 
80  PRINT 
90  PRINT 
100  PRINT 
110  INPUT 
120  COLOR 
130  CLS 
140  CLS 
150  COLOR 
1 


,  1 

"THE  QUATECH  BOARD  IS  A  12  BIT  A/D  CONVERTER” 

"IT  IS  SETUP  FOR  +  AND  -  2.5  VOLTS  INPUT  MAX" 

"THE  SAMPLE  RATE  AND  NUMBER  OF  DATA  POINTS" 

"ARE  USER  SELECTABLE" 

"PRESS  RETURN  TO  CONTINUE  " ,CCC$ 

,0 

,1 

6 


0 


'*********************************************************** 


170  'THE  FIRST  TEN  LINES  OF  THE  PROGRAM  ARE  FOR  LINKING  THE 
BASIC  * 

180  'CALL  COMMANDS  TO  THE  QUATECH  MACHINE  LANGUAGE 


1  9  0 

************************************************************ 

200  ADC . SETUP=&H3 : SETCTM=&H6 
SETC0=&H9 

I NADC1 2 . B=&HC : SEGADDR=&HF : I NADC1 2 . BB=&H 1 2 
SETMD 1 =&H 1 5 : SETMD2=&H 18:SETC1=&H1B:SETC2=&H1E 
RSETC1 =&H2 1 : RSETC2=&H24 : READC1 =&H27 : READC2=&H2A 
DEF  SEG=0 

CSEG2=PEEK ( &H3CA )+256*PEEK( &H3CB ) 

DEF  SEG=CSEG2 

GOSUB  330  'QUATECH  BOARD  SETUP 

GOSUB  860  'DATA  FILE  SETUP 

GOSUB  1010  ' AQU I  RE  DATA 


210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 


END 

************************************************** 
THIS  IS  THE  SUBPROGRAM  THAT  SETS  UP  THE  QUATECH  * 
BOARD  FOR  SAMPLE  RATE  AND  NUMBER  OF  DATA  POINTS.* 
************************************************** 


DIM  ARRAY%(9000) 

DIM  VOLT! (9000) 
ADDRESS%=&H300 
N%  =  2 

CALL  ADC . SETUP ( ADDRESS^ , N% ) 


'BOARD  ADDRESS 
'#  OF  CHANNELS 
SETS  BOARD  ADD 


PR  I  NT 
PR  I  NT 


'CHOSE  A  SAMPLE  RATE  BETWEEN  200  HZ  -  32000  HZ" 
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470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 


INPUT  "  INPUT  SAMPLE  RATE  IN  Hz:  ",HZ* 

CO*  =  5000000! /HZ*  *  SETS  SAMPLE  RATE 

CALL  SETCO  (CO*) 

PR  I  NT 
PR  I  NT 

PR  I  NT 

1 

PRINT  "CHOSE  THE  NUMBER  OF  DATA  POINTS  BETWEEN  65-8191” 
PR  I  NT 

INPUT  "ENTER  NUMBER  OF  DATA  POINTS:  "  ,  LENGTH* 

COLOR  ,0 
CLS 

COLOR  ,  1 


PRINT  "SAMPLE  RATE:  "  HZ* 

t 

PR  I  NT 
PR  I  NT 
PR  I  NT 

l 

PRINT  "NUMBER  OF  DATA  POINTS:  "  LENGTH* 

t 

PR  I  NT 
PRINT 
PRINT 

TM  =  ( 1 /HZ* ) *LENGTH* 

PRINT  TM  "  SEC  TOTAL  DATA  ACQUISITION  TIME" 

PR  I  NT 
PR  I  NT 
PR  I  NT 

INPUT  "ARE  THE  SAMPLE  RATE  AND  NUMBER  OF  POINTS  CORRECT 
V/N  ";T$ 

IF  T$=  "N"  OR  T$  =  "n"  THEN  470 
» 


COLOR  ,0 
COLOR  , 1 
RETURN 

' ****************************************************** 

'THIS  SUBPROGRAM  WILL  SETUP  YOUR  DATA  FILE.  * 

'****************************************************** 

COLOR  ,0 
CLS 

COLOR  ,  1 

t 

INPUT  "NAME  DATA  FILE  FOR  DATA  STORAGE:  XXXXX  NO 
EXTENSION,  5  CHR  OR  LESS  ";NA$ 

« 

NAT$=  "D:"  +  NA$  +  "T"  +  ".DAT" 

' NA$  =  "D:"  +  NA$  +  "F"  +  ".DAT” 
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980  'OPEN  NA$  FOR  OUTPUT  AS  #1 
990  OPEN  NAT$  FOR  OUTPUT  AS  #2 
1000  RETURN 

10 10  '*************************************************** 
1020  ’  THIS  SUBPROGRAM  IS  TO  ACQUIRE  DATA  AND  CHANGE  IT  * 
1030  '  TO  A  VOLTAGE.  * 

1040  ' *************************************************** 


"PRESS  RETURN  TO  ACQUIRE  DATA. 


1010  ’******************************** 

1020  ’  THIS  SUBPROGRAM  IS  TO  ACQUIRE  D 
1030  '  TO  A  VOLTAGE. 

1040  ' ******************************** 
1050  ' 

1060  COLOR  ,0 
1070  CLS 
1080  COLOR  , 1 

1090  INPUT  "PRESS  RETURN  TO  ACQUIRE  DA 
1100  ’ 

1110  COLOR  ,0 
1 1 20  CLS 
1 130  COLOR  , 1 

1140  PRINT  "ACQUIRING  DATA  !!!!!!!!!!" 
1150  TM  =  ( 1/HZ%)*LENGTH% 

1155  PRINT 

1160  PRINT  TM  "  SEC  TOTAL  DATA  ACQU I 
1170  ' 

1180  CALL  I NADC 1 2 . B (LENGTHS , ARRAYS ( 1 ) ) 

1184  PRINT 

1185  PRINT  "DATA  HAS  BEEN  ACQUIRED  AND 


COLOR  ,0 
CLS 

COLOR  , 1 

PR  I  NT  "ACQU I R I NG  DATA  !!!!!!!!!!” 

TM  =  ( 1/HZ%)*LENGTH% 

PR  I  NT 

PRINT  TM  "  SEC  TOTAL  DATA  ACQUISITION 


1185 
1  190 
1200 
1210 
1220 
1230 
1235 
1240 
1250 
1260 
1265 
1270 
1280 
1290 
1300 
1310 
1320 


PRINT  "DATA  HAS  BEEN  ACQUIRED  AND  I S  BE  I 
T!  =0 
TT  =  0 

FOR  I  =  1  TO  (LENGTHS*2)  STEP  2 
T I =  1/HZ%  *  TT 

VOLT ! ( I ) =  - ( ARRAYS ( I ) + . 5 ) *  (5/4095)+  2.5 
VOLT • ( I + 1 ) =  - ( ARRAYS ( I + 1 ) + . 5 ) *  (5/4095)+ 
' PR  I NT#1 , USING"##. ####";  VOLT! ( I ) 
PRINT#2, USING"##. ###~~~*  T !  ; 
PRINT#2,USING"##.###  ";  VOLT ! ( I ) ; 

PR  I NT#2 , US  I NG"## . ###" ;  VOLT ! ( I  +  1 ) 

TT=TT+1 
NEXT  I 
COLOR  ,0 
CLS 

COLOR  , 1 
RETURN 


’STORE  INTO  ARRAY 


BEING  STORED. 


This  program  was  written  in  Microsoft  GW-Basic  version  3.1 
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MATLAB  PSD  and  Digital  Filtering  Program 


echo  on 
c  1  c 

%  This  program  filters  data  through  a  Chebyshev  lowpass 
%  and  stopband  filter.  Then  the  PSD  is  plotted  after  each 
%  filtering. 

pause  %  Strike  any  key  to  continue, 
c  1  c 

n= 1 ength ( y  )  ; 

plot(y( 1 :n) ) ,  title('lnput  Data'),  pause 
c  1  c 

N  =  4;  %  Filter  order 

stopband  =  [.0225  .0255];  %  Stopband  specification 

[Bb,Ab]  =  butter(4,  stopband ,' stop ') ; 

stopband  =  [.023  .025];  %  Stopband  specification 

ripple  =  40;  %  Allowable  attenuation  in  decibels 

[Be, Ac]  =  cheby2(N,  ripple,  stopband stop ') ; 

N  =  8;  %  Filter  order 

passband  =  [.95];  %  Lowpass  specification 

[Bbl.Abl]  =  butter(N,  passband); 

ripple  =  .01;  %  Allowable  ripple,  in  decibels 

[Bel, Acl]  =  cheby1(N,  ripple,  passband); 

%  Let's  filter  the  data  with  the  Chebyshev  stopband  filter 
%  and  plot  the  PSDF . 

Y  =  f i ltf i lt(Bc,Ac,y)  ; 

P  =  spectr um( Y , 1 024 ) ; 

speep lot(P , 5000) ,  grid,  pause 
c  1  c 

%  Let's  filter  the  data  again  with  the  Chebyshev  lowpass 

%  filter  and  plot  the  resulting  PSD. 

Z  =  f i ltf i lt(Bc1 , Acl , Y) ; 

Q  =  spectrum( Z , 1 024 ) ; 

specplot (Q, 5000) ,  grid,  pause 
c  1  c 


This  program  was  written  using  MATLAB  version  3.5f. 
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Sample  LOTUS  123  Cell  Formulas  for  Reducing  Reset  Data 


A  2 :  0.0008 
B 2  :  0.564 
C2:  -2.501 
02:  ( A  1  - $A$ 1 ) 

E2:  (§1 F(@ABS(C2)>2 ,@TRUE , ©FALSE )  =  1  ) 
F2:  (®IF(E2=1 ,  0 , 0 . 95 1 *C2 ) ) 

G2:  ( F2 ) 

H2:  (G2-12. 492*02) 

A3:  0.001 
B  3 :  0.612 
C3:  0.153 
03:  ( A2-$A$ 1 ) 

E3:  ( @ I F ( @ABS (C3 ) >  2 , ©TRUE , @FALSE ) = 1 ) 
F  3 :  (@IF(E3=1 ,0,0.951*03)) 

G3:  (F3) 

H3:  (G3- 1 2 . 492*D3 ) 

A4:  0.0012 
B4 :  -0.23 
C4 :  0.171 
04:  ( A3-$A$ 1 ) 

E4:  ( @ I F (@ABS ( C4 ) >  2 , @TRUE , ©FALSE )  =  1  ) 
F  4 :  (@iF(E4=1 ,0,0.951 *C4 ) ) 

G4:  ( F4 ) 

H4 :  (G4-12. 492*04) 

A5:  0.0014 
B5 :  -2.121 
C5:  0.079 
D5 :  ( A4-$A$ 1 ) 

E5:  (§IF(@ABS(C5)>2,@TRUE,@FALSE)  =  1  ) 
F5:  (§IF(E5=1 ,0,0.951*05)) 

G5:  (F5) 

H5:  (05-12.492*05) 

A6:  0.0016 
B6 :  -0.882 
C6:  0.211 
06:  ( A5-$A$ 1 ) 

E6:  ( ® I F (@ABS (C6 ) >  2 , ®TRUE ,@FALSE)  =  1 ) 

F  6 :  (@IF(E6=1 ,0,0. 951 *C 6)) 

G6:  (F6) 

H6 :  (Q6-12. 492*D6) 

A7:  0.0018 
B7:  -1.86 
C7:  0.079 
07:  ( A6-$A$ 1 ) 

E7:  (® I F (@ABS (C7 ) >  2 , ©TRUE , ©FALSE ) = 1 ) 

F  7 :  (®IF(E7=1 ,0,0. 951*C7)) 

07:  (F7) 

H7 :  (07-12.492*07) 
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A8 :  0.002 
B8 :  -0.464 
C8:  0.087 
D8 :  (A7-$A$1) 

E8:  (@IF(@ABS(C8)>2,§TRUE,@FALSE)  =  1  ) 

F  8 :  (@IF(E8=1 , 0 , 0 . 95  1  *C8  )  ) 

G8 :  @ I F(E8=0 , (F8+$G$7 ) , ( ( $G$7 +G7 +F 9 ) / 2 ) ) 

H8 :  (G8-12.492*D8) 

A9 :  0.0022 
B9:  0.294 
C9 :  0.105 
D9:  (A8-$A$1) 

E  9 :  (@l F (@ABS(C9) >2, ©TRUE, ©FALSE )=1 ) 

F  9 :  (@IF(E9=1 ,0,0.951*09)) 

G9:  @IF(E9=0, (F9+$G$7) , ( ( $G$7 +G8 +F 1 0 ) / 2 ) ) 

H9:  (G9-12.492*D9) 

A  1  0 :  0.0024 
BIO:  -0.674 
CIO  :  0.122 
DIO:  (A9-$A$1) 

E 1 0 :  (©IF (@ABS ( Cl  0 ) >  2  , ©TRUE , @F ALSE ) = 1 ) 

F 1 0 :  (@IF(E10=1 , 0 , 0 . 95 1 *C1 0 ) ) 

G 1 0 :  ©  I  F ( E  1  0  =  0 , (F10+$G$7) , (($G$7+G9+F1 1 )/2)) 

H 1 0 :  (G10-12.492*D10) 

All:  0.0026 
B  1  1  :  -1.724 
C 1  1  :  0.018 
Dll:  (A10-$A$1 ) 

Ell:  (@ I F (@ABS (C 1 1 ) >  2 , @TRUE , @F  ALSE )  =  1 ) 

F  1 1  :  (§IF(E1 1  =  1 ,0,0.951*C1 1 )) 

ail:  @IF(E11=0,(F1 1 +$G$7 ) , ( ( $G$7 +G 1 0+F 1 2 ) / 2 ) ) 

H 1 1 :  (G1 1-12. 492*D1 1 ) 

A 1  2  :  0.0028 
B 1 2 :  -1.895 
C 1 2 :  -2.501 
D 1 2 :  (All -$A$ 1 ) 

El  2:  (@ I F ( @ABS ( C 1 2 ) >  2 , ©TRUE , @FALSE )  =  1 ) 

F 1 2 :  (@IF(E12=1 ,0,0.951 *C 12)) 

G 1 2 :  §IF(E12=0, (F12+$G$1 1 ) , (($G$1 1+G1 1+F13)/2)) 
H 1 2 :  (G12-12.492*D12) 

A 1 3 :  0.003 
B 1 3 :  -2.303 
C 1 3 :  0.029 


D 1  3 : 
E  1  3 : 
F  1  3 : 
G  1  3 : 
H 1  3 : 


( A 1 2-$A$ 1 ) 

(@!F(@ABS(C13)>2,@TRUE,@FALSE)=1) 
(@IF(E13=1 ,0,0. 951 *C 13)) 

@IF(E13=0, (F13+$G$1 1 ) , (($G$1 1 +G 1 2+F 1 4 ) /2 ) ) 
(G1 3- 12.492*013) 


A  1 4 :  0.0032 


B 1 4 :  -0.326 


C 1 4 :  0.017 


D 1 4 :  ( A  1 3-$A$ 1 ) 

E 1 4 :  (@l F(@ABS(C14)>2 , §TRUE , ©FALSE )= 1 ) 
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F 1 4 :  (@IF(E14=1 ,0,0.951 *C 14)) 

G 1 4 :  @ I F ( E 1 4-0 , (F14+$G$1 1 ) , ( ($G$1 1 +G 1 3+F 1 5 ) / 2 ) ) 
H 1 4 :  (G14-12.492*D14) 

A  1 5 :  0.0034 
B 1 5 :  -0.248 
C 1 5 :  0.033 
D 1 5 :  ( A  1 4-$A$ 1 ) 

E 1 5 :  (@IF(@ABS(C15)>2,@TRUE,@FALSE)  =  1  ) 

F  1  5 :  (@IF(E15=1 ,0,0.951*015)) 

G 1 5 :  @ I F ( E 1 5=0 , ( F 1 5+$G$ 1 1 ) , ( ($G$1 1+G14+F16)/2)) 
H 1 5 :  (G15-12.492*D15) 

A  1 6 :  0.0036 
B 1  6 :  0.913 
C 1 6 :  -0.023 
D 1 6 :  (A  1 5-$A$1 ) 

E 1 6 :  (@ I F ( @ABS ( C 1  6  )  >  2  , @TRUE , ©FALSE )  =  1 ) 

F 1 6 :  (@IF(E16=1 ,0,0.951*016)) 

G 1 6 :  @IF(E16=0, (F16+$G$1 1 ) , ( ( $G$ 1 1 +G 1 5+F 1 7 ) / 2 ) ) 
H 1 6 :  (G16- 1 2 . 492*D1 6 ) 

A  1 7 :  0.0038 
B 1  7 :  -0.614 
C 1  7  :  0.021 
D 1  7  :  ( A 1 6-$A$1 ) 

El  7:  (@ I F (@ABS(C1 7 ) >2 , ©TRUE , ©FALSE )  =  1  ) 

F 1 7 :  (@IF(E17=1 ,0,0.951*017)) 

G17:  @IF(E17  =  0, (F17  +  $G$ 1 1 ) , ( ($G$1 1 +G1 6  +  F 1 8 )/2 ) ) 
H 1 7 :  (G17-12. 492*017) 

A  1 8 :  0.004 
B 1 8 :  0.371 
C  1  8 :  0.009 
D  1 8 :  ( A  1 7-$A$ 1 ) 

E 1 8 :  (@IF(@ABS(C18)>2,@TRUE,@FALSE)=1 ) 

F 1 8 :  (§IF(E18=1 ,0,0.951*018)) 

G 1 8 :  @IF(E18=0, (F18+$G$1 1 ) , ( ($G$1 1+G1 7+F 1 9 )/2) ) 
H 1 8 :  (G18-12.492*D18) 

A 1  9 :  0.0042 
B  1  9 :  1.943 
C 1  9 :  -0.07 
D 1 9 :  ( A 1 8-$A$ 1 ) 

El  9:  (@ I F ( @ABS ( C 1 9 ) >  2 , @TRUE , ©FALSE )  =  1 ) 

F 1 9 :  (§IF(E19=1, 0,0. 951*019)) 

G 1 9 :  @!F(E19=0 , (F19+$G$1 1 ) , ( ($G$ 1 1 +G 1 8+F20 ) /2 ) ) 
H 1 9 :  (G1 9- 1 2 . 492*D 1 9 ) 

A20:  0.0044 
B20:  1.748 
C20 :  0.011 
020 :  ( A 1 9-$A$ 1 ) 

E  20 :  (@l F(§ABS(C20)>2 , ©TRUE , @FALSE )  =  1 ) 

F20 :  (§IF(E20=1 ,0,0.951*020)) 

G20 :  @!F(E20=0, (F20+$G$1 1 ) , ( ( $G$ 1 1 +G 1 9+F2 1 ) /2 ) ) 
H20 :  (G20-12.492*D20) 
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A  2 1 :  0.0046 
B21:  2.499 
C21:  -0.029 
D  2 1 :  (A20-$A$1) 

E  2 1 :  (@ I F ( @ABS ( C2 1 ) >  2 , @TRUE ,@FALSE)=1) 

F  2 1 :  (@IF(E21=1, 0,0. 951*021)) 

G21 :  @1 F(E21=0, (F21+$G$1 1 ) , ( ( $G$ 1 1 +G20+F 2 2 ) / 2 ) ) 
H21:  (G2 1  -  1 2 . 492*D2  1  ) 

A  22 :  0.0048 
B22:  0.245 
C22:  -2.501 
D22:  ( A2 1  - $A$ 1  ) 

E22 :  (@IF(@ABS(C22)>2,@TRUE,@FALSE)  =  1  ) 

F  2  2 :  (@IF(E22=1 ,0,0.951*022)) 

G22:  §1 F(E22=0 , (F22+$G$21 ) , ( ( $G$2 1 +G2 1 +F 2 3 ) / 2 ) ) 
H22:  (G22-12.492*D22) 

A23:  0.005 
B23:  -1.339 
C23 :  -0.076 
D23:  (A22-$A$1) 

E23 :  (@ I F (@ABS ( C2 3 ) >  2 , @TRUE , ©FALSE )  =  1 ) 

F23:  (@IF(E23=1 ,0,0.951*023)) 

G23:  @IF(E23  =  0, (F23+$G$21 ) , ( ( $G$2 1 +G2  2  +  F24 ) /2 ) ) 
H23 :  (G23-12.492*D23) 
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Table  D-1.  PVDF  Film  Transducer  Cost 


Mater  i  a  1 _ _ 

PVDF  Film 

AD549  Operational  Amplifier 
3N170  Transistor 
Twin  Lead  Shielded  Cable 
AD525  Instrumentation  Amplifier 
8  Channel  Multiplexing  Chip 


Cost _ 

$  0.20  cm* 

$  17.00  each 
$  3.00  each 

$  0.20  ft 

$  25.  00  each 


0 . 42  each 


Approximate  cost  for  eight  transducers: 


$  60.00 


Table  D-2.  Dynamic  Pressure  Transducer  Cost 


Mater i a  1 


Cost 


Endevco  8510-B  Pressure  Transducer 


$  350.00  each 


Signal  Conditioner 


350.00  each 


Approximate  cost  for  eight  transducers 


$5600.00 
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Appendix  E:  Digital  Filter  Data 


MATLAB  Filter  Demonstration  Program 


echo  on 
c  1  c 

%  This  example  shows  the  use  of  the  FFT  function  for 
%  spectral  analysis.  A  common  use  of  FFT's  is  to  find  the 
%  frequency  components  of  a  signal  buried  in  a  noisy  time 
%  domain  signal . 


pause  %  Strike  any  key  to  continue, 
c  1  c 

%  First  we  need  to  create  some  data.  Consider  data  sampled 
%  at  5000  Hz.  We  start  by  forming  a  time  axis  for  our 
%  data,  running  from  t  =  0  until  t=.41  in  steps  of  0.2 
%  m i 1 1 i second : 

t  =  0: .0002: .4096; 

pause  %  Strike  any  key  to  continue. 

%  Next  we  can  form  a  signal  containing  60  Hz: 
x  =  s i n( 2*p i *60*t )  ; 

pause  %  Strike  any  key  to  continue. 

%  and  add  some  random  noise  with  a  standard  deviation  of  2 
%  to  produce  a  noisy  signal  y: 

rand( ’ normal ' ) 
y  =  x  +  2*rand(t); 

pause  %  Strike  any  key  to  continue. 

X=[ ] ; t=[ ] ; 

c  1  c 

%  Let's  take  a  look  at  our  noisy  signal  y  by  plotting  it. 
pause  %  Strike  any  key  for  plot. 

plot(y( 1 : 1024) ) ,  title('Noisy  time  domain  signal'),  pause 
c  1  c 

%  Clearly,  it  is  difficult  to  identify  the  frequency 

%  components  from  looking  at  the  original  signal;  that's 

%  why  spectral  analysis  is  so  popular. 

%  Finding  the  discrete  Fourier  transform  of  the  noisy 

%  signal  y  is  easy.  We  just  take  the  fast-Fourier  transform 

%  (FFT);  using  the  Spectrum  Function: 
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P  =  spectrum(y , 1024) ; 

pause  %  Strike  any  key  to  continue. 

%  The  power  spectral  density,  a  measurement  of  the  energy 
%  at  various  frequencies,  is  found  with: 

%  P =spectr um( y , 1 024 ) 

pause  %  Strike  any  key  to  continue. 

c  1  c 

%  To  plot  the  power  spectral  density,  we  must  first  form  a 
%  frequency  axis: 

f  =  5000/2048* ( 0  :  1024)  ; 

%  which  we  do  for  the  first  1024  points.  (The  remainder  of 
%  the  1024  points  are  symmetric.)  We  can  now  plot  the 
%  power  spectral  density: 

pause  %  Strike  any  key  for  plot. 

specp 1 ot (P , 5000 ) ,  pause 
c  1  c 

%  Let's  zoom  in  and  plot  only  up  to  600  Hz: 
pause  %  Strike  any  key  for  plot. 

plot(f(1:128),P(1:128,1)),  t i t 1 e( ’ Power  spectral  density'),  .. 
x 1 abe 1 (' Frequency  (Hz)'),  pause 

rand( ' uni  form' ) 
c  1  c 

%  This  example  demonstrates  the  use  of  the  BUTTER  and 
%  CHEBY2  functions  to  design  stopband  filters.  A  sampling 

%  rate  of  5000  Hz  is  assumed. 

%  We  start  by  specifying  the  desired  frequency  response, 

%  point-wise,  with  1.0  correspond i ng  to  half  the  sample 

%  rate: 

f  =  [0  .0236  .0236  .0244  .0244  1]; 

H  =  [  1  1  0  0  1  1  ]  ; 

pause  %  Strike  any  key  to  continue. 

c  1  c 

%  Next  we  plot  the  desired  frequency  response  to  make  sure 

%  it  is  what  we  want.  To  plot  the  frequency  response,  we 

%  can  scale  the  normalized  frequency  axis: 
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fs  =  5000;  %  sampling  rate 

fhz  =  f*fs/2; 

%  Now  we  are  ready  to  plot, 
pause  %  Strike  any  key  for  plot. 

plot(fhz.H),  t i t 1 e( ' Des i red  Frequency  Response’),., 
x 1 abel (' Frequency  (Hz)’),  y 1 abe 1 (' Magn i tude ') ,  pause 
c  1  c 

%  If  the  desired  frequency  response  did  not  look  right, 

%  we  would  go  back  and  change  H.  However,  in  this  case 
%  everything  looks  OK  so  we  perform  the  design: 

format  compact 

n  =  1024; 

ff  =  fs/(2*n)  *  ( 0 : n- 1  )  ; 

%  Now  let’s  design  Butterworth  and  Chebyshev  stopband 

%  filters 

%  with  the  same  stopband  (defined  between  0.0  and  1.0). 

N  =  4  ;  'v  r  i  1  ter  order 

stopband  =  [.0225  .0255];  %  Stopband  specification 

[Bb,Ab]  =  butter(4,  stopband ,’ stop ’)  ; 

stopband  =  [.023  .025];  %  Stopband  specification 

ripple  =  40;  %  A’lowable  attenuation  in  decibels 

[Be, Ac]  =  cheby2(N,  ripple,  stopband ,’ stop ’)  ; 

%  Now  find  the  frequency  responses: 

hb  =  f r eqz ( Bb , Ab , n )  ; 
he  =  f r eqz ( Be , Ac , n ) ; 
hm  =  [ abs ( hb ) ]  ; 
hn  =  [ abs ( he  )  ]  ; 

pause  %  Strike  any  key  for  plot. 

plot(ff,hm),  t it  1 e( * Butterworth  filter’),  .. 
x 1 abe 1 (’ Frequency  (Hz)’),  y 1 abe 1 (’ Magn i tude ’) ,  grid,  pause 
c  1  c 

%  Often  it  is  useful  to  look  at  the  frequency  response  on 
%  logarithmic  (dB)  scales. 

pause  %  Strike  any  key  for  plot. 

plot(ff( 15: 35) , 20*1 og10(hm( 15 : 35 ,:))),. . 
t it  1 e( ’ Butterworth  filter’),  .. 

x 1 abe 1 (’ Frequency  (Hz)’),  y 1 abe 1 (’ Magn i tude  in  dB’),grid, 
pause 
c  1  c 

plot(ff,hn),  t i t 1 e( ’ Chebyshev  filters'),  .. 
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x 1 abe 1 (' Frequency  (Hz)'),  y 1 abe 1 ( ' Magn i tude ' ) ,  grid,  pause 
c  1  c 

%  Often  it  is  useful  to  look  at  the  frequency  response  on 
%  logarithmic  (dB)  scales. 

pause  %  Strike  any  key  for  plot. 

plot(ff(  15: 35) ,20*log10(hn( 15:35, :))),.  . 
t i t 1 e ( ' Cheby shev  filters'),  .. 

x 1 abe 1 (' Frequency  (Hz)'),  y 1 abe 1 ( ' Magn i tude  in  dB'),grid, 
pause 
c  1  c 

%  This  example  demonstrates  the  use  of  the  BUTTER  and 
%  CHEBY1  functions  to  design  low-pass  filters.  A  sampling 

%  rate  of  5000  Hz  is  assumed. 

%  We  start  by  specifying  the  desired  frequency  response, 
%  point-wise,  with  1.0  cor r espond i ng  to  half  the  sample 
%  rate: 

f  =  [0  .95  .95  1]; 

H  =  [  1  1  0  0]  ; 


c  1  c 

%  Next  we  plot  the  desired  frequency  response  to  make  sure 
%  it  is  what  we  want.  To  plot  the  frequency  response,  we 

%  can  scale  the  normalized  frequency  axis: 

%  Now  we  are  ready  to  plot. 

fhx  =  f *f s/2 

pause  %  Strike  any  key  for  plot. 

plot(fhx.H),  t i t 1 e( ' Des i red  Frequency  Response'),., 
x 1 abe 1 (' Frequency  (Hz)'),  y 1 abe 1 (' Magn i tude ’) ,  grid,  pause 
c  1  c 

%  If  the  desired  frequency  response  did  not  look  right, 

%  we  would  go  back  and  change  H.  However,  in  this  case 
%  everything  looks  OK  so  we  perform  the  design: 

format  compact 

c  1  c 


%  Now  let's  design  Butterworth  and  Chebyshev  low-pass 
%  filters  with  the  same  band  (defined  between  0.0  and  1.0). 


N  =  8; 
passband 
[ Bb 1 , Ab 1 ] 


%  Filter  order 

[.95];  %  Low-pass  specification 

butter(N,  passband); 


ripple  =  .01;  %  Allowable  ripple,  in  decibels 

[Bel, Acl]  =  cheby1(N,  ripple,  passband); 
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%  Now  find  the  frequency  responses: 

hbl  =  freqz(Bb1 ,Ab1 ,n)  ; 
hcl  =  freqz(Bc1 ,Ac1 ,n) ; 
hi  =  [abs(hbl)  abs(hcl)]; 

pause  %  Strike  any  key  for  plot. 

p 1 ot ( f f , h 1  )  ,  t i 1 1 e( ' Butterworth  and  Chebyshev  filters’),  .. 
x 1 abe 1 (' Frequency  (Hz)'),  y 1 abe 1 (' Magn i tude '), pause 
c  1  c 

%  Often  it  is  useful  to  look  at  the  frequency  response  on 
%  logarithmic  (dB)  scales. 

pause  %  Strike  any  key  for  plot. 

p lot (ff(2:n), 20* log  1 0 ( h  1  ( 2  :  n  ,  :))),.  . 

t i t 1 e( ’ Butterworth  and  Chebyshev  filters'),  .. 

x 1 abe 1 (' Frequency  (Hz)'),  y 1 abe 1 ( ' Magn i tude  in  dB ’  )  ,  grid, 

pause 

c  1  c 

%  Let's  filter  the  data  with  the  Chebyshev  stopband  filter 
%  and  plot  the  PSD. 

Y  =  f i ltf i lt(Bc,Ac,y) ; 

P  =  spectrum(Y, 1024) ; 

specplot(P , 5000) ,  grid,  pause 

c  1  c 

%  Let's  filter  the  data  again  with  Chebyshev  low-pass 

%  filter  and  plot  the  resulting  PSD. 

Z  =  f i 1 tf i 1 1 ( Be  1 , Ac  1 , Y  )  ; 

Q  =  spectrum( Z , 1 024 )  ; 

specplot(Q,5000) ,  grid,  pause 
c  1  c 
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i  E-1.  Random  Noise  Signal  with  Embedded  60  Hz  Sine  Wave 


Figure  E-2.  PSD  of  the  Unfiltered  Random  Noise  Signal. 
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Figure  E-4.  Response  of  the  Butterworth  Stopband  Filter. 
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Figure  E-6.  Response  of  the  Butterworth  and  Chebyshev  Low-pass 
F i 1 ters . 
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Figure  E-7.  PSD  of  the  'Noisy'  Data  Filtered  with  the 
Chebyshev  Stopband  Filter. 


Figure  E-8.  PSD  of  the  'Noisy'  Data  Filtered  the  Chebyshev 
Stopband  and  Low-pass  Filters. 
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This  appendix  gives  a  sample  of  how  the  data  was 
reconstructed  for  the  sine  wave  shown  in  Figure  4-6.  Before 
any  calculations  are  done,  the  offset,  due  to  an  improper 
nulling  adjustment  on  the  charge  amplifier,  is  removed.  This 
is  done  by  subtracting  the  simple  mean  of  all  the  data  points. 
For  point  [1:1],  the  adjusted  value  is  0.0055  -  0.0003  = 
0.0053,  where  0.0003  is  the  mean  offset.  The  first  three 
adjusted  data  sample  periods  are  shown  in  the  following  arrays 
(see  Table  F-1  column  E).  That  is, 

Array  1  =  [0.0055  0.0085  0.0135  0.0205  0.0275  0.0365  0.0475 
0.0595  0.0755  0.0885  0.1075] 

Array  2  =  [0.0145  0.0325  0.0545  0.0765  0.0985  0.1205  0.1475 
0.1715  0.1945  0.2235  0.2495] 

Array  3  =  [0.0155  0.0425  0.0725  0.0975  0.1295  0.1555  0.1865 
0.2145  0.2415  0.2685  0.2955]. 

To  create  time  history,  the  first  point  of  Array  1  [1:1] 
is  taken  as  the  starting  point  and  set  to  zero.  The  rest  of 
the  sample  points  remain  unaffected. 

The  first  point  in  Array  2  [2:1],  is  a  "reset"  point  and 
is  discarded.  It  is  replaced  by  a  "bridge"  point  calculated 
from  equation  3-36.  That  is, 

[2:1]  -  xn-i+xn-i+xn*i  -  0 . 1023  * 0 . 1023  ♦  0 . 0309  -0  1177 
2  2 

where  x,  is  the  point  [2:1],  x,.(  is  point  [1:11],  and  x,,(is 
point  [2:2].  The  next  ten  points  in  Array  2,  now  need  to  have 
the  value  of  point  [1:11]  added  to  them.  That  is, 


[2:2],,,  = 

[2:2],  „  + 

[1:11] 

[2:3],,,  = 

[  2  :  3],l4  + 

[1:11] 

[2:11],.,  = 

[2:H], „  + 

[1:11] 

The  "new"  points  replace  the  "old"  values.  This  procedure  is 
repeated  for  the  next  array  of  points  (see  Table  F-1 
column  I).  The  resulting  new  points  are  shown  in  the 
following  arrays.  That  is, 

Array  1  =  [0.0053  0.0081  0.0129  0.0195  0.0262  0.0347  0.0452 
0.0566  0.0718  0.0842  0.1023] 

Array  2,„  =  [0.1177  0.1332  0.1541  0.1750  0.1960  0.2169  0.2426 
0.2654  0.2875  0.3148  0.3396] 

Array  3,„  =  [0.3597  0.3800  0.4085  0.4323  0.4627  0.4874  0.5170 
0.5436  0.5693  0.5949  0.6206]. 

Next,  the  DC  offset  is  removed  by  scaling  all  points  with  the 
DC  offset  factor.  For  example  point  [1:2]  becomes  0.95  * 
0.0081  =  0.0077.  All  the  other  points  are  scaled  in  the  same 
manner  (see  Table  F-1  column  H). 

In  the  case  of  the  sine  wave,  the  rise  rate  was  found  by 
fitting  a  linear  least-squares  line  through  the  peak  amplitude 
points  of  each  cycle.  These  three  values  (0.916,  0.896,  and 
0.868)  were  0.0196  seconds  apart.  The  least-squares  analysis 
gave  a  slope  of  -1.173. 

The  rise  rate  error  was  removed  by  subtracting  the  rise  or 
"slope  correction  factor"  from  each  point.  Point  [1:1]  was 
the  starting  point  and  is  assumed  to  have  no  error.  Point 
[1:2]  was  corrected  in  the  following  manner: 
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=  M:2]  +  1/5555  *  n  *  1.173  =  0.0079 
where  n  is  the  index  beginning  with  n  =  0  for  [1:1],  n  =  12 
for  [2:1],  n  =  22  for  [3:1],  and  so  forth  for  the  remaining 
points  in  the  data  record.  Finally,  the  sine  wave  was 
re-adjusted  by  subtracting  0.4646  from  each  point  to  give  a 
zero  mean  to  the  corrected  data  points  (see  Table  F-1  column 
J). 

For  the  first  thirty  three  points,  the  final  corrected 
values  are  as  follows:  [-0.4646  -0.4563  -0.4513  -0.4446 

-0.4376  -0.4288  -0.4181  -0.4066  -0.3911  -0.3785  -0.3602 

-0.3446  -0.3289  -0.3077  -0.2866  -0.2655  -0.2443  -0.2185 

-0.1954  -0.1733  -0.1456  -0.1206  -0.1002  -0.0797  -0.0510 

-0.0270  0.0036  0.0286  0.0583  0.0851  0.1110  0.1369  0.1628]. 

These  values  correspond  to  the  first  thirty  three  points 
plotted  in  Figure  4-6. 
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Table  F-1.  Sample  Numbers  Used  to  Reconstruct  Sine  Data. 
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0 . 00036 
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-0 
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-0 
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-0 
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'  0 

-  0 .034739 

0.034739 

-0.42882 

0.00108 

0.047529 

-0 

65928 

1  0 

'  0 . 0452 

:  0.0452 

-0.41314 

0.00126 

0 .059529 

-0 

.73828 

;  o 

j  0.056612 
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-0.40652  : 
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-0 

81228 

0 
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0 
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1 

0 
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0.002161 
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0 
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0.133196 
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0 . 002341 
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15028 

0 

0.051857 

0.1541  18 

-0.30775 
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-1 

20728 

0 

0.072779 

0 . 17504 

-0.28661 

0.002701 
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-1 
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0 

0.093701 

0. 195962 

-0.26548  | 

0 . 00288 1 

0. 120529 

- 1 

31028 

0 

0 .  1  14623 

0. 216884 

-0.24435  ! 

0 . 00306 1 

0.147529 

- 1 

35528 

0 

0.1403 

0 . 242561 

-0.21346  . 

0 . 00324 1 

0. 171529 

- 1 

39428 

0 

0 . 163124 

0 . 265385 

-0.  19542 

0 . 00342 1 

0 . 194529 

- 1 

43128 

0 

0 . 184997 

0 . 287258 

-0. 17334  j 

0.003601 

0 . 223529 

- 1 
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0 
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0.003781 

0 . 249529 
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0 
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1 

0 
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- 1 

52828 

0 

0 . 040445 

0.380008 
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0.004321 

0.07252° 

1 

54128 

0 
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0.408538 
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0 . 00450  1 

0.097529 

- 1 
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0 

0.09275 
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0.004681 
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0 

0. 123182 

0 . 462745 

0.003624 
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0. 155529 

- 1 

55028 

0 

0 . 147908 

0.487471 

0 . 02856 1 

0 . 005041 

0. 186529 

- 1 
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0 
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0 . 005221 
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- 1 
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0.204017 

0.54358 

0.085092 

0.005401 

0.241529 

-1 

51328 

0 

0.229694 

0.569257 

0. 110981 

0 . 00558 1 

0.268529 

-1 
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0 

0.255371 

0.594934 

0 . 136869 

0.005721 
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- 1 . 

46228 

0 

0.281048 

0.620611 

0 . 162757 

0.005941 

0.011529 

-1 

43228 

1 

0 

0 . 639407 
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Column  A 
Co  1 umn  B 
Column  C 
Co  1 umn  D 
Co  1 umn  E 
Column  F 
Column  G 
Column  H 
Co  1 umn  I 
Co  1 umn  J 


Not  shown  -  Original  time  tag. 

Not  shown  -  Original  reset  data. 

Not  shown  -  Original  input  data. 

Time  adjusted  to  zero. 

Reset  data  offset  adjusted  to  zero. 

Input  data  offset  adjusted  to  zero. 

Reset  point  logic:  0=not  reset,  1=reset. 

DC  offset  removed  from  reset  data. 

Reset  data  joined  to  form  time  history. 

Reset  data  corrected  for  rise  rate  and  zero  mean. 
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LOTUS  123  Cell  Formulas  Used  to  Reconstruct  the  Sine  Data 


A  1  :  0.00324 
B  1  :  -0.01 
Cl  :  -0.14 
D 1 :  ( A  1  - $A$ 1 ) 

El:  ( B 1 -$B$343 ) 

FI:  (C1-$C$343) 

G  1  :  0 

HI:  (@IF(G1=1 ,0,0. 951*E1 )) 

11:0 

J 1 :  (1. 1  7267  303*01 +  ( 11-0.46461)) 

A2 :  0.00342 
B2:  -0.007 
C2 :  -0.227 
D2:  (A2-$A$1) 

E  2 :  (B2-$B$343) 

F2 :  (C2-$C$343) 

G 2:  (@IF(@ABS(E2-E1 ) >0.051 , @TRUE , ©FALSE )= 1 ) 
H2 :  (@IF(G2=1 ,0,0.951*E2)) 

12:  §IF(G2=0,(H2*1) , ( 11+1 1+H3)/2)) 

J2:  (1 . 17267303*D2+( 12-0.46461)) 

A3:  0.0036 
B3 :  -0.002 
C3 :  -0.314 
D3:  (A3-$A$1) 

E  3 :  ( B3-$B$343 ) 

F  3 :  (C3-$C$343) 

G3 :  (@ I F (@ABS (E3-E2 ) >0 . 05 1 , ©TRUE , ©FALSE )= 1 ) 
H3 :  (@IF(G3=1 ,0,0.951*E3)) 

I  3:  @1 F(G3=0, (H3*1 ) , ( I  2+ I 2+H4)/2) ) 

J3:  (  1 . 17267303*D3+( I  3-0. 46461  )) 

A4 :  0.00378 
B4 :  0.005 
C4  :  -0.394 
D4:  ( A4-$A$1 ) 

E4 :  ( B4-$B$343) 

F4 :  (C4-$C$343 ) 

G4 :  (@ I F (@ABS( E4-E3 ) >0.051 , ©TRUE , ©FALSE )= 1 ) 
H4 :  (@IF(G4=1 ,0,0 . 951*E4) ) 

14:  @IF(G4=0, (H4*1 ) , ( i 3+1 3+H5)/2)) 

J4 :  (1 . 17267303*D4+( 14-0.46461)) 

A5 :  0.00396 
85:  0.012 
C5 :  -0.479 
D5 :  (A5-$A$1) 

E5 :  ( B5-$B$343 ) 

F  5 :  (C5-$C$343) 

G5:  (©IF (@ABS ( E5-E4 ) >0 . 05 1 , ©TRUE , @FALSE ) = 1 ) 
H5:  (§IF(G5=1 , 0 , 0 . 95 1 *E5  )  ) 

15:  ©IF(G5=0,(H5*1) , ( I  4+ I 4+H6 ) /2 ) ) 

J5:  (  1 . 17267303*D5  +  ( I  5-0. 46461  )) 


F-5 


A6 :  0.00414 
B6 :  0.021 
C6 :  -0.559 
D6 :  (A6-$A$1) 

E  6 :  (B6-$B$34  3) 

F  6 :  (C6-$C$343 ) 

G6 :  (@IF(@ABS(E6-E5)>0.051 , ©TRUE , @F ALSE )  =  1  ) 
H6:  (@!F(G6=1,0,0.951*E6)) 

16:  © I F (G6  =  0 , (H6* 1 ) , ( I  5+ I 5  +  H7  )/2  )  ) 

J6 :  (1 . 17267303*D6+( 16-0.46461)) 

A7 :  0.00432 
B7 :  0.032 
C7 :  -0.641 
D7 :  ( A7 -$A$ 1 ) 

E7 :  (B7-$B$343) 

F7 :  (C7-$C$343) 

G7 :  (@l F(@ABS(E7-E6)>0 . 051 , ©TRUE , ©FALSE )= 1 ) 
H7:  (@IF(G7=1 , 0 , 0 . 95 1 *E7  )  ) 

17:  @IF(G7  =  0, (H7*1  )  ,  ( I  6+ I 6  +  H8  )  / 2 ) ) 

J7 :  (  1 . 17267303*D7  +  ( I  7-0 . 46461  )  ) 

A8:  0.0045 
B8:  0.044 
C8 :  -0.72 
08 :  (A8-$A$1 ) 

E8 :  ( B8-$B$343 ) 

F8 :  (C8-$C$343) 

G8 :  (@ I F(@ABS(E8-E7 ) >0 . 05 1 , ©TRUE , ©FALSE )= 1 ) 

H8 :  (@IF(G8=1 ,0,0.951*E8)) 

18:  @IF(G8=0,(H8*1),( I  7+ 1 7  +H9 ) / 2 ) ) 

J8:  (  1 . 17267303*D8  +  ( I  8-0.46461  ) ) 

A9 :  0.00468 
B9 :  0.06 
C9 :  -0.794 
D9 :  (A9-$A$1) 

E9 :  ( B9-$B$343 ) 

F  9 :  (C9-$C$343) 

G9 :  (@IF(@ABS(E9-E8)>0.051 , ©TRUE , ©FALSE )= 1 ) 

H9 :  (©IF(G9=1,0,0.951*E9)) 

19:  @IF(G9=0, (H9*1 ) ,(!8+l8+Hl0)/2)) 

J9:  (1 . 17267303*D9+( 19-0.46461)) 

A 10:  0.00486 
BIO:  0.073 
CIO:  -0.867 
DIO:  ( A  1 0-$A$ 1 ) 

E 1 0 :  ( B 10-$B$343 ) 

F 1 0 :  (C10-$C$343) 

G 1 0 :  (@IF(@ABS(E10-E9)>0.051 , ©TRUE , ©FALSE )= 1 ) 
H 1 0 :  (@IF(G10=1 ,0,0.951*E10)) 

110:  «IF(G10=0,(H10*1), (l9+l9+H11)/2)) 

J 1 0 :  (  1 . 1 7267303*D10+( I  10-0.46461  )) 

All:  0.005041 
B 1 1 :  0.092 
Cl  1 :  -0.939 
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D 1  I 
Ell 
FI  1 
G 1  1 
HI  1 
I  1  1 
J1  1 
A  1  2 
B  1  2 
C 1  2 
D  1  2 
El  2 
F  1  2 
G 1  2 
H  1  2 
I  12 
J1  2 
A  1  3 
B  1  3 
Cl  3 
013 
E  1  3 
F  1  3 
G 1  3 
H  1  3 
I  13 
J  1  3 
A  1  4 
B  1  4 
C 1  4 
D  1  4 
E  1  4 
F  1  4 
G 1 4 
H  1  4 
I  14 
J  1  4 
A  1  5 
B  1  5 
Cl  5 
D 1  5 
E  1  5 
F  1  5 
G15 
H 1  5 
I  15 
J 1  5 
A 1  6 
B  1  6 
C 1  6 
016 
E  1 6 
F  1  6 


(All -$A$ 1 ) 

( B 1 1 -$B$34  3 ) 

(Cl  1 -$C$343 ) 

(@IF(@ABS(E1 1-E 10) >0.051 , ©TRUE , ©FALSE )  =  1  ) 
(@ I F (G 1 1  =  1 , 0 , 0 . 95 1 *E 1  1  )) 

@IF(G11=0,(H11*1), (110+1  1 0  +  H 1 2 ) /2  )  ) 

(.  1  .  17  267303*01  1  +  (  I  1  1-0.46461  )) 

0.005221 
-0.001 
-1.007 
( A 1 2-$A$ 1 ) 

( B 1 2-$B$343 ) 

(Cl  2- *C$343) 

(@l F(@ABS(E 12-El  1 ) >0.051 , ©TRUE , ©FALSE )  =  1 ) 
(@l F(G1 2= 1 , 0 , 0 . 95 1 *E 1 2  ) ) 

@IF(G12=0,(Hl2+$l$11),(($l$11+l11+H13)/2)) 
(1.17267 30 3*0 1 2  +  ( I  12-0. 46461)) 

0 . 005401 
0.017 
-1 .073 
( A  1 3-$A$ 1 ) 

( B 1 3-$B$34  3 ) 

(Cl  3- $C$343 ) 

(@IF(@ABS(E13-E12)>0.051 , ©TRUE , ©FALSE )= 1 ) 
(§IF(G13=1 ,0,0.951*E13) ) 

©I F(G13  =  0, (H 1 3  +  $ I  $  1 1 ) , (($1 $1 1+ I  1 2  +  H 1 4 ) /2 ) ) 
(  1 .  1 7  267303*01 3+ ( I  13-0.46461  )) 

0.005581 
0.039 
-1.132 
( A  1 4-$A$ 1 ) 

( B 1 4-$B$343 ) 

(C14-$C$343) 

(©IF (@ABS (E14-E13)>0.051 , ©TRUE , ©FALSE )= 1 ) 
(@ I F (G 1 4= 1 ,0,0.951*E14)) 

©I F(G14  =  0, (H14  +  $l $1 1 ) , ( ($l$1 1+ I  1 3+H 1 5 ) / 2 ) ) 
(  1 . 17267303*D14+( I  14-0.46461  ) ) 

0.005761 
0.061 
-1.189 
(A  1 5-$A$ 1  ) 

(B15-$B$343) 

(C15-$C$343) 

(@IF(@ABS(E15-E14)>0.051 , ©TRUE , ©FALSE )= 1 ) 
(@IF(G15=1 , 0 , 0 . 95 1 *E  1  5  ) ) 

©I F(G15=0, (H15+$l $1 1 ) , ( ($1 $1 1+ I 14+H16)/2) ) 
(  1 . 1 7  267303*0 1 5  + ( I  15-0.46461  )) 

0.005941 
0.083 
-1.243 
( A  1 6-$A$ 1 ) 

( B 1 6-$B$343 ) 

(Cl 6-$C$343 ) 
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G 1 6 :  (@IF(@ABS(E 16-El  5) >0.051 , ©TRUE ,@FALSE)=1) 

H 1 6 :  (@IF(G16=1 ,0,0.951*E16)) 

116:  @IF(G16=0,(H16+$l$11),(($l$11+l15+H17)/2)) 
J 1 6 :  (  1 .  17267303*D16  +  ( I  16-0.46461  )) 

A  1  7 :  0.006121 
B 1  7 :  0.105 
C 1 7 :  -1.292 
D 1 7 :  ( A  1 7  - $A$ 1 ) 

E 1 7 :  (B17-$B$343) 

F 1 7 :  (Cl  7 -$C$343 ) 

G 1 7 :  (@IF(@ABS(E 17-El  6) >0.051 , ©TRUE , ©FALSE )  =  1 ) 

H 1 7 :  (@IF(G17=1 ,0,0.951*E17)) 

I  17:  @IF(G17=0,(H17+$I$11),(($I$11+I  16  +  H18)/2) ) 
J 1 7 :  ( 1 . 1 7267303*01 7+ ( 117-0.46461)) 

A  1 8 :  0.006301 
B  1  8 :  0.132 
C 1 8 :  -1.337 
018:  (A  1 8-$A$ 1  ) 

E 1 8 :  ( B 1 8-$B$343 ) 

F 1 8 :  (Cl 8-$C$343 ) 

G 1 8 :  (@IF(@ABS(E18-E17)>0.051 , ©TRUE , ©FALSE )= 1 ) 

H 1 8 :  (@IF(G18=1  , 0 , 0 . 95 1 *E 1 8  )  ) 

118:  @IF(G18=0,(H18+$l$11),(($l$11+l17+H19)/2)) 
J 1 8 :  (  1 . 17  26  7303*01 8  + ( I  18-0.46461)) 

A  1 9 :  0.006481 
B 1 9 :  0.156 
C 1 9 :  -1.376 
019:  ( A  1 9-$A$ 1 ) 

E  1  9 :  ( B 1 9-$B$343 ) 

F 1 9 :  (C 1 9-$C$343 ) 

G1 9 :  (@IF(@ABS(E19-E 18) >0.051 , ©TRUE , ©FALSE )  =  1  ) 

H 1 9 :  (@IF(G19=1 ,0,0.951*E19)  ) 

119:  @IF(G19=0,(H19+$I$11),(($I$11+I 18+H20)/2)) 
J 1 9 :  (  1 . 17  267303*01 9  + ( I  19-0.46461  ) ) 

A20:  0.006661 
B20 :  0.179 
C20 :  -1.413 
020:  ( A20-$A$ 1 ) 

E20 :  (B20-$B$343) 

F20 :  (C20-$C$343) 

G20:  (©IF(@ABS(E20-E19)>0.051 , ©TRUE , ©FALSE )= 1 ) 
H20 :  (@IF(G20=1 ,0,0. 951*E20) ) 

120:  @IF(G20=0,(H20+$I$11),(($I$11+I 19+H21)/2)) 
J20:  (  1 . 1 7267303*D20  +  ( I  20-0.46461  ) ) 

A21:  0.006841 
B  2 1 :  0.208 
C21:  -1.444 
021:  ( A2 1 -$A$ 1 ) 

E21:  ( B2 1 -$B$343 ) 

F  2 1 :  (C2 1 -$C$343) 

G21 :  (©IF(@ABS(E21-E20)>0.051 ,@TRUE ,@FALSE)=1) 
H21:  (@IF(G21=1 ,0,0.951*E21)) 

121:  @IF(G21=0,(H21+$I$11),(($I$11+I 20+H22)/2)) 
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J21 :  (  1 . 17  267  30  3*D2 1  +  ( 121-0.46461  )) 

A22:  0.007021 
B22 :  0.234 
C22:  -1.474 
D22:  (A22-$A$1) 

E  2  2 :  (B22-$B$343) 

F  2  2 :  (C22-$C$343) 

G22:  (@ I F(@ABS(E22-E21 )> 0.051 , ©TRUE , ©FALSE )= 1 ) 
H22:  (@IF(G22=1 ,0,0.951*E22)) 

122:  @IF(G22  =  0, (H22  +  $!$1 1 ) , (($1 $1 1  +  I 21+H23)/2)  ) 
J22:  (  1 . 17  267  303*D22+ ( 122-0.46461  )) 

A23:  0.007201 
B23:  0 
C23 :  -1.496 
D23:  (A23-$A$1) 

E  2  3 :  (B23-$B$343) 

F23:  (C23-$C$343) 

G23:  (@ I F(@ABS(E23-E22 ) >0 . 05 1 , ©TRUE , ©FALSE )= 1 ) 
H23:  (@IF(G23=1 , 0 , 0 . 95  1  *E23 ) ) 

I  23:  @1 F (G23=0 , (H23+$ I $22) , ( ( $ t  $2  2+ I 22  +  H24 ) /2 ) ) 
J2  3:  (  1 . 1 7  267303*023+ ( I  23-0.46461  )) 

A24:  0.007381 
B24:  0.027 
C24 :  -1.51 
D24 :  (A24-$A$1) 

E24 :  ( B24-$B$343 ) 

F24 :  (C24-$C$343) 

G24 :  (@  I F (@ABS ( E24-E23 ) >0.051 , ©TRUE , ©FALSE )= 1 ) 
H24 :  (©IF (G24= 1 ,0,0.951*E24)) 

I  24:  @IF(G24  =  0, (H24  +  $l$22) , ( ( $  I $2  2+ I 23+H25)/2)) 
J24:  (  1 . 1 7  267  303*024+ ( I  24-0.4646  I  )  ) 

A25 :  0.007561 
B25:  0.057 
C25 :  -1.523 
D25 :  (A25-$A$1) 

E25 :  ( B25-$B$343 ) 

F25 :  (C25-$C$343 ) 

G25 :  (@IF(@ABS(E25-E24)>0.051 ,@TRUE , ©FALSE )= 1 ) 
H25 :  (@IF(G25=1 , 0 , 0 . 95 1 *E25  )  ) 

I  25:  @IF(G25  =  0, (H25  +  $l$22) , ( ($1 $22+ I 24  +  H26)/2) ) 
J25 :  (  1 . 1 7  26  7  303*D25  +  ( I  25-0.46461  )) 

A26 :  0.007741 
B26:  0.082 
C26 :  -1.532 
026:  (A26-$A$1 ) 

E26 :  (B26-$B$343 ) 

F  26 :  (C26-$C$343) 

G26 :  (@IF(@ABS(E26-E25)>0.051 , ©TRUE , ©FALSE )= 1 ) 
H26 :  (@IF(G26=1 ,0,0.951*E26)) 

I  26:  @IF(G26  =  0, (H26  +  $l$22) , ( ($1 $22+ I 25+H27)/2)) 
J26:  (  1 . 1 7267303*D26  +  ( I  26-0.46461  ) ) 

A27 :  0.007921 
B27 :  0.114 
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C27 :  -1.54 
D27 :  ( A2  7 -$A$ 1 ) 

E  2  7 :  (B27-$B$343) 

F27 :  (C27-$C$343) 

G2  7:  (@IF(@ABS(E27-E26)>0.051 , ©TRUE , §FALSE ) = 1 ) 

H27 :  (@IF(G27=1 , 0 , 0 . 95 1 *E27  )  ) 

127:  @IF(G27=0, (H27+$l$22) , ( ($i$22+ I 26+H28)/2) ) 
J27:  (  1 . 17  267 30  3*D27  +  ( I  27-0.46461  )) 

A28 :  0.008101 
B28 :  0.14 
C28 :  -1.532 
D28 :  (A28-$A$1) 

E28 :  (B28-$B$343) 

F  28 :  (C28-$C$343) 

G28 :  (@IF(@ABS(E28-E27)>0.051 , @TRUE , ©FALSE )= 1 ) 
H28:  (@IF(G28=1 ,0,0.951*E28)) 

I  28:  @IF(G28=0,(H28+$l$22),(($l$22+l27+H29)/2)) 
J28:  (  1 . 17267 303*D28+ ( I  28-0.46461  )) 

A29:  0.008281 
B29:  0.171 
C29 :  -1.521 
D29:  (A29-$A$1) 

E  29 :  ( B29-$B$343 ) 

F29 :  (C29-$C$343 ) 

G29 :  (@IF(@ABS(E29-E28)>0.051 , ©TRUE , ©FALSE )  =  1  ) 
H29 :  (@IF(G29=1 , 0 , 0 . 95 1 *E29 ) ) 

I  29:  @IF(G29  =  0, (H29  +  $l$22) , ( ( $  I $22+ I 28+H30 ) /2 ) ) 
J29:  (1 . 1 7267303*029+ ( 129-0.46461)) 

A30 :  0.008461 
B30 :  0.199 
C30 :  -1.509 
D30 :  ( A30-$A$ 1 ) 

E  30 :  ( B30 -$B$343 ) 

F  30 :  (C30-$C$343 ) 

G30 :  (@IF(@ABS(E30-E29)>0.051 , ©TRUE , ©FALSE )  =  1  ) 
H30:  (§IF(G30=1 , 0 , 0 . 95 1 *E 30 ) ) 

I  30:  @IF(G30  =  0,(H30  +  $l$22),(($l$22+l29+H31)/2)) 
J30 :  (  1 . 1 7267303*D30  +  ( I  30-0.46461  ) ) 

A3 1 :  0.008641 
B31:  0.226 
C31:  -1.495 
D31:  ( A3 1 -$A$ 1 ) 

E31:  (B31-$B$343) 

F31:  (C31-$C$343) 

G31 :  (@IF(@ABS(E31-E30)>0.051 , @TRUE , §FALSE )  =  1  ) 
H31:  (§IF(G31=1 , 0 , 0 . 95 1 *E3 1 ) ) 

131:  §IF(G31=0,(H31+$l$22),(($t$22+l30+H32)/2)) 
J31 :  (  1 . 1 7  267  303*D3 1  +  ( 131-0.46461  )) 

A32:  0.008821 
B32:  0.253 
C32 :  -1.474 
D32 :  (A32-$A$1) 

E  32 :  ( B  32 -$B$343 ) 
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F  32 :  (C32-$C$343) 

G32:  (§IF(@ABS(E32-E31 )> 0.051 , ©TRUE , ©FALSE )= 1 ) 

H32:  (@IF(G32=1 , 0 , 0 . 95 1 *E32  )  ) 

I  32:  @1 F(G32  =  0, (H32  +  $l $22) , ( ($1 $22+ I 31+H33)/2) ) 
J32 :  ( 1 . 1 7267303*032+ ( 132-0.46461)) 

A  33 :  0.009001 
B  33 :  0.28 
C33 :  -1.444 
033:  (A33-$A$1) 

E33 :  (B33-$B$343) 
m  •  ( c  ■sa 3 1 

G33 :  (@ I F (@ABS ( E  33-E32 ) >  0 . 05 1 , ©TRUE , ©FALSE )  =  1 ) 

H33:  (@IF(G33=1 , 0 , 0 . 95 1 *E33) ) 

I  33:  ©I F(G33=0, (H33+$l $22) , ( ($1 $22+ I 32  +  H34)/2) ) 
J33:  (  1 . 17 26 7 303*033+ ( I  33-0.46461  )) 

A34:  0.009181 
B34 :  -0.004 
C34 :  -1.414 
D34 :  ( A34-$A$ 1 ) 

E34 :  (B34-$B$343) 

F  34 :  (C34-$C$343) 

G34:  (@IF(@ABS(E34-E33)>0.051 , ©TRUE , ©FALSE )= 1 ) 

H34:  (@IF(G34=1 ,0,0.951*E34)) 

134:  ©  I F (G34  =  0 ,(H34  +  $l$33),(($l$33+l33  +  H35)/2)) 
J34:  (  1  .  17 26 730 3*034+ ( 134-0.46461  )) 
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The  purpose  of  this  research  was  to  investigate  the 
char acter i st i cs  of  a  small  and  inexpensive  pressure  transducer  made  from 
piezoelectric  film.  The  pressure  transducer  was  designed  to  be  small  and 
surface  mountable  so  that  the  pressure  distribution  could  be  measured  on 
an  aerodynamic  body.  A  charge  reset  technique  was  used  to  overcome  the 
traditional  limitations,  low  frequency  measurements  and  charge  decay,  of 
a  piezoelectric  device.  The  resulting  time  history  had  to  be  integrated, 
so  a  method  was  developed  to  integrate  the  data  and  remove  sources  of 
error.  The  time  history  reconstr uct i on  technique  was  tested  with  a  known 
input  and  found  to  accurately  reassemble  the  original  input.  The 
piezoelectric  film  pressure  transducer  was  tested  in  a  flow  field  along 
side  a  commercial  pressure  transducer.  The  piezoelectric  transducer’s 
frequency  response  showed  agreement  with  the  frequency  response  of  the 
commercial  transducer.  Also,  the  time  history  response  of  the 
piezoelectric  transducer  gave  evidence  that  low  frequency  pressure 
measurements  are  possible  with  a  piezoelectric  pressure  transducer. 

Pressure  Transducer,  Pressure  Gage,  Piezoelectric  Film,  143 
PVDF ,  Aerodynamic  Testing,  Pressure  Measurement 
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